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7-10 GALAXIES IN THE HUDF AND GOODS FIELDS, AND THEIR UV LUMINOSITY FUNCTIONS 
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ABSTRACT 

We use all available deep optical ACS and near-IR data over both the HUDF and the two GOODS 
fields to search for star-forming galaxies at z > 7 and constrain the UV LF within the first 700 Myrs. 
Our data set includes ~ 23 arcmin 2 of deep NICMOS J + H data and ~ 248 arcmin 2 of ground-based 
(ISAAC+MOIRCS) data, coincident with ACS optical data of greater or equal depths. In total, 
we find 8 <z>^7.3 z-dropouts in our search fields, but no z ~ 9 ./-dropout candidates. A careful 
consideration of a wide variety of different contaminants suggest an overall contamination level of just 
~ 12% for our z-dropout selection. After performing detailed simulations to accurately estimate the 
selection volumes, we derive constraints on the UV LFs at z ~ 7 and z ~ 9. For a faint-end slope 



+ 1.7 



x 10" 



a = —1.74, our most likely values for M^ v and <f>* at z ~ 7 are —19.8 ± 0.4 mag and 1.1 
Mpc~ 3 , respectively. Our search results for z ~ 9 J-dropouts set a la lower limit on M^ v of —19.6 
mag assuming that 0* and a are the same as their values at slightly later times. This lower limit on 
My V is 1.4 mag fainter than our best-fit value at z ~ 4, suggesting that the UV LF has undergone 
substantial evolution over this time period. No evolution is ruled out at 99% confidence from z ~ 7 
to z ~ 6 and at 80% confidence from z ~ 9 to z ~ 7. The inferred brightening in M^y with redshift 
(i.e., M* v = (-21.02 ± 0.09) + (0.36 ± 0.08) (z- 3.8)) matches the evolution expected in the halo mass 
function, if the mass-to-light ratio of halos evolves as ~ (1 + z)~ x . Finally, we consider the shape 
of the UV LF at z > 5 and discuss the implications of the Schechter-like form of the observed LFs, 
particularly the unexpected abrupt cut-off at the bright end. 
Subject headings: galaxies: evolution — galaxies: high-redshift 



1. INTRODUCTION 

One of the most important goals of observational cos- 
mology has been to identify galaxies at the earliest times 
and to characterize their evolution. Over the past few 
years, substantial progress has been made on this front, 
particularly at high redshifts (z ^4 — 6). Much of this 
progress has been due to the significant quantity of deep 
multi- wavelength data made possible by instruments like 
the Hubble Space Telescope (HST) Advanced Camera 
for Surveys (Ford et al. 2003) or Subaru's Suprime-Cam 
(Miyazaki et al. 2002). As a result, current selections 
(e.g., Yoshida et al. 2006; Beckwith et al. 2006; Bouwens 
et al. 2007) now number in the tens of thousands of galax- 
ies and reach to luminosities as faint as —16 AB mag 
(~0.01 L* z=3 ). This has enabled astronomers to measure 
a wide-range of properties in high-redshift galaxies and 
examine their evolution across cosmic time. 

One of the key observables considered in previous stud- 
ies of high-redshift galaxies is the luminosity function in 
the rest-frame UV. The UV LF is of significant interest 
not only because of the quantitative constraints it im- 
poses on early galaxy formation, but also because of the 
information it provides on the star formation rate and 
UV output of galaxies at the earliest times. This latter 
information is important for investigating the role galax- 
ies may have played in the reionization of the universe 



(e.g., Stanway et al. 2003; Giavalisco et al. 2004b; Yan 
& Windhorst 2004; Bouwens et al. 2006). 

In principle, the large statistical samples of star- 
forming galaxies at z ~ 4 — 6 just described allow us 
to determine the UV LFs and their evolution in great 
detail out to z ~ 6. Unfortunately, despite the large 
sample sizes and wide luminosity ranges of current selec- 
tions, there continues to be some uncertainty regarding 
how the UV LF evolves at very high redshift (i.e., from 
z > 6 to z ~ 3). Some groups (Shimasaku et al. 2005; 
Bouwens et al. 2006; Yoshida et al. 2006; Ouchi et al. 
2004a; Bouwens et al. 2007; Dickinson et al. 2004) have 
argued the evolution occurs primarily at bright end of 
the LF, others (Iwata et al. 2003; Sawicki & Thompson 
2006; Iwata et al. 2007) have asserted it happens more at 
the faint-end, and still others have suggested the evolu- 
tion is independent of luminosity (Beckwith et al. 2006). 
The diversity of conclusions drawn by different teams 
would seem to suggest that systematics are likely playing 
a large role. A careful discussion of this issue as well as 
detailed comparisons with contemporary determinations 
of the UV LF at z ~ 4 — 6 can be found in Bouwens et 
al. (2007). 

One way of gaining additional leverage in looking at 
the evolution of the UV LF at high redshift is to push 
to even higher redshift, i.e., z ~ 7 — 10, to see what 
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galaxies were like at these epoches. Since it is expected 
that changes in the UV LF will be much more substantial 
over the interval z ~ 7 ~ 10 to z ^ 3 — 4 than from 
z ~ 6 to z ~ 4, it should be relatively easy to ascertain 
how the UV LF is evolving at early cosmic times. The 
challenge, of course, is that the number of star-forming 
galaxy candidates at z ~ 7 — 10 has been very small, of 
order ~ 4 — 5 galaxies (Bouwens et al. 2004c; Bouwens 
et al. 2005; Bouwens & Illingworth 2006; lye et al. 2006; 
Labbe et al. 2006; Mannucci et al. 2007; Stark et al. 
2007b; Bradley et al. 2008). As a result, the bulk of 
the effort has been directed towards simply establishing 
the overall prevalence of bright star-forming galaxies at 
z ~ 7 - 10 (Bouwens et al. 2004c; Yan & Windhorst 
2004; Bouwens et al. 2005; Bouwens & Illingworth 2006; 
Mannucci et al. 2007; lye et al. 2006; Stanway et al. 
2008). 

In the present work, we will try to take studies of the 
UV LF at z > 7 one step further and not just look at 
the overall prevalence of star-forming galaxies at z > 7, 
but at how the volume density of these galaxies appears 
to vary as a function of luminosity. In doing so, we will 
attempt to derive constraints on the overall shape of the 
UV LF at z > 7 and therefore evaluate different mod- 
els for how the UV LF evolves with cosmic time. In 
particular, we are interested in looking at how the vol- 
ume density of the bright population evolves and how 
this compares to the evolution of the faint population. If 
the mean luminosity of galaxies really does increase with 
time at early cosmic times, as one might expect from 
hierarchical models, the volume density of bright galax- 
ies should evolve much more rapidly at z > 7 than the 
volume density of the faint population. This change has 
been alternatively parametrized as a brightening of the 
characteristic luminosity M^ v with cosmic time (e.g., 
Bouwens et al. 2006; Yoshida et al. 2006; Bouwens et al. 
2007) or a flattening of the faint-end slope a (e.g., Yan 
& Windhorst 2004). 

To help with this endeavor, we will make use of the 
previous data used by our group in high-redshift (z ~ 
7—10) galaxy searches (e.g., Bouwens & Illingworth 2006; 
Bouwens et al. 2005) and also take advantage of a sig- 
nificant amount of newer data. Among these new data 
are deep ACS observations over the NICMOS parallels to 
the HUDF (Oesch et al. 2007), a large number of moder- 
ately deep NICMOS observations over the two GOODS 
fields (H. Teplitz et al. 2008, in prep), the deep wide- 
area ISAAC observations over the CDF-South GOODS 
field (Mannucci et al. 2007; Stanway et al. 2008; B. Van- 
dame et al. 2008, in preparation; J. Retzlaff et al. 2008, 
in prep), and the deep wide-area MOIRCS observations 
over the HDF-North GOODS field (Kajisawa et al. 2006; 
Ouchi et al. 2007). We will also be taking advantage of 
improvements in our NICMOS reduction techniques to 
push our previous selections slightly fainter. These new 
data (particularly the new wide-area ground-based data) 
will be crucial for examining the UV LF at z > 7 over a 
wider range in luminosity than was possible in previous 
work (e.g., Bouwens et al. 2004c; Bouwens & Illingworth 
2006; Mannucci et al. 2007; Stanway et al. 2008), pro- 
viding us with additional leverage to constrain the shape 
of the UV LF at > 7. The ultra deep ACS data over 
the NICMOS parallels to the HUDF will also allow us to 
increase both the depth and robustness of our searches 



for z ~ 8 — 10 galaxies. 

We will organize this paper as follows. In §2, we will 
provide a brief description of the observations we will use 
in this paper. In §3, we will describe how we select our 
dropout samples from the imaging data and how we es- 
timate the likely contamination levels. In §4, we will use 
these samples to derive the rest-frame UV LF at z ~ 7 
and set constraints on this LF at z ~ 9. Finally, in §5, we 
discuss our results, and in §6 we will provide a summary 
of the main conclusions from this work and provide a 
short outlook for future progress. Throughout this work, 
we will find it convenient to denote luminosities in terms 
of the characteristic luminosity £* =3 at z ~ 3 (Steidel 
et al. 1999), i.e., Muv,ab = —21.07. Where necessary, 
we assume f2 = 0.3, fl\ = 0.7, H = 70km/s/Mpc. Al- 
though these parameters are slightly different from those 
determined from the WMAP five-year results (Dunkley 
et al. 2008), they allow for convenient comparison with 
other recent results expressed in a similar manner. Our z 
and J dropout criteria select galaxies over a range in red- 
shift (i.e., z ~ 6.5 — 8.0 and z ~ 8 — 10, respectively), with 
mean redshifts <z>=7.3 and <z>~ 9.0, respectively (see 
§4.2). For succintness, we will usually refer to the red- 
shift ranges selected by these dropout criteria as z <~ 7 
and z <~ 9. We will express all magnitudes in the AB 
system (Oke & Gunn 1983). 

2. OBSERVATIONAL DATA 

A summary of the observational data used to search 
for z > 7 z and J dropouts is provided in Table 1. 

2.1. NICMOS Search Fields with Deep zJH coverage 

Our deepest, wide-area constraint on the rest-frame 
UV LF at z ~ 7—10 comes from the numerous HST 
NICMOS NIC3 fields which lie in and around the two 
GOODS fields. Our reductions of the NICMOS data 
over and around the GOODS fields (including the HUDF 
NICMOS data) was performed with the python package 
"nicred.py" (Magee et al. 2007) and represent a slight 
update to the reductions we performed for the Bouwens 
& Illingworth (2006) work. The most important change 
in our reduction procedure is that we now include a cor- 
rection for the non-linear response of the NIC3 camera 
(de Jong et al. 2006). This has the effect of brightening 
the measured magnitudes in the Jn and Hi 60 bands by 
~ 0.15 and ~ 0.05 mag (relative to what these magni- 
tudes would be without the corrections). Our reductions 
of the optical ACS data are the same as those considered 
in our most recent work on the UV LF at z ~ 4 — 6 
(Bouwens et al. 2007) and include essentially all primary 
and parallel ACS data taken in the vicinity of the two 
GOODS fields (e.g., Giavalisco et al. 2004a; Bouwens et 
al. 2004a; Riess et al. 2007; Oesch et al. 2007). 

Our primary NICMOS data set includes a somewhat 
larger number of search fields, some of significant depth, 
and better optical coverage than was used in our previous 
search for z ~ 7 galaxies (Bouwens & Illingworth 2006: 
see Figure 1 from the Supplementary Information) . The 
most notable additions are (1) one ultra deep (28.1 AB 
mag at 5a) ~ 0.7 arcmin 2 NIC3 pointing within the area 
of the HUDF (Oesch et al. 2007), (2) ultra deep (28.7 
AB mag at 5a) optical ACS coverage over the second 
NICMOS parallel to the HUDF (Oesch et al. 2007), and 
(3) <~ 4 arcmin 2 of additional deep (26.7 AB mag at 
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Fig. 1. — Deep near-IR fields in and around the CDF-South (left) and HDF-N (right) GOODS fields where we conducted searches for 
star-forming galaxies at z > 7 using z and J dropout criteria. Areas with deep (> 28 mag at 5<r) or ultra-deep (> 29 mag at 5<r) optical 
V606*775 2 850 coverage are shaded in blue and dark blue, respectively. Areas with deep (> 26.5 mag at 5<r) or ultradeep (> 28 mag at 5<x) 
NICMOS coverage are shown in orange and red, respectively. For clarity, we have labelled each of these search fields with a letter, "U" 
indicating the HUDF Thompson field (Thompson et al. 2005), "S" indicating the HUDF Stiavelli field (Oesch et al. 2007), "D" indicating 
the HDF-North Dickinson field (Dickinson 1999), "T" indicating the HDF-North Thompson field, (Thompson et al. 1999), "1" indicating 
HUDF-NICPAR1 (Bouwens et al. 2005; Oesch et al. 2007), "2" indicating HUDF-NICPAR2 (Bouwens et al. 2005; Oesch et al. 2007), and 
"G" representing the GOODS parallel fields (S. Malhotra et al. 2008, in prep; H. Teplitz et al. 2008, in prep). The deep, wide-area VLT 
ISAAC coverage of the CDF-South GOODS field is shown as a light orange region. The ~ 112 arcmin 2 of deep MOIRCS coverage (the 
4 GTO pointings) over the HDF-North GOODS field is shown with a similar orange shading (the deepest MOIRCS pointing having the 
darkest shading). Additional information on our search fields can be found in Table 1. 



5a) NICMOS data within the GOODS area (Riess et al. 
2007; H. Teplitz et al. 2008, in preparation). 

We have also made improvements to our reductions 
of the two deep NICMOS fields taken in parallel to 
the HUDF (Bouwens & Illingworth 2006; Oesch et al. 
2007: see Figure 1). Reducing the NICMOS data over 
these fields is challenging because of the small number 
of unique dither positions used in the acquisition of the 
data and the large number of individual frames (~ 200 
frames in each field). As a result, any systematics present 
in characterizing the properties of individual pixels on 
the NIC3 detector propagate directly into the reduced 
frames, and therefore it becomes necessary to calibrate 
these properties sufficiently well so that systematics do 
not have a dramatic effect on the overall noise charac- 
teristics. To accomplish this, we generated super median 
images from a stack of all the NIC3 frames taken over 
a field, subtracted this median image from each NIC3 
frame, and then redrizzled the individual frames together 
to produce a cleaner reduction. Before generating the 
super median images, we mask out all sources identi- 
fied in the x 2 image (Szalay et al. 1999) created from 
our initial reductions of the Jno and Hiqq band imaging 
data; ~ 300 — 400 sources were masked per NIC3 frame. 
We elected to use the x 2 image for this process to en- 
sure that the same sources were masked in both the Jno 
and Hiqo frames. Therefore, the colors of sources were 
not affected by this "median stacking and subtraction" 
procedure. To obtain the best possible reduction of the 
data, we repeated the above procedure several times: us- 
ing each successive reduction of the data to better mask 
out the faint sources in the data, median stacking the 
masked frames, subtracting the medians from the indi- 
vidual frames, and then drizzling the data together to 



generate a cleaner "final" reduction. 

In the dropout searches we performed over the NIC- 
MOS data, we only included data which satisfied certain 
minimal depth requirements in both the optical and near- 
IR data to avoid the edges of the individual reductions 
where the number of individual frames contributing to 
any given pixel was small (i.e., where there might be 
concerns about some non-Gaussian characteristics in the 
noise). In our z ~ 7 zgso-dropout search, we only in- 
cluded areas where the zsso, >/iio> and -ffi6o depths were 
27.5 AB mag, 26.8 AB mag, and 26.5 AB mag, respec- 
tively, over a 0.6"-diameter aperture (5a). We demanded 
greater depth for our Jno-band imaging data than we 
did for the -Hi6o-band data given the importance of a 
robust detection in the Jno-band to making the case for 
an object being at z > 7. For our z ~ 9 Jno-dropout 
search, we only included areas where the zsso, Jno, and 
ffi6o-band depths were 27.5 AB mag, 26.8 AB mag, and 
26.8 AB mag, respectively, over a 0.6"-diameter aperture 
(5a). We set more stringent requirements on the Hiqq- 
band depths for our Jno dropout search because of the 
greater need to use the i/i6o-band data to establish the 
reality of the sources. In total, the area of our search 
fields were ~ 23 arcmin 2 for our zsso-dropout selection 
and ~ 21 arcmin 2 for our Jno-dropout selection. 

2.2. ISAAC Search Fields 

The deep (J ~ 25.5 AB mag, H,K ~ 25.0 AB mag 
at 5a), wide-area (~ 143 arcmin 2 ) ISAAC data over the 
CDF-South GOODS field (B. Vandame et al. 2008, in 
prep; J. Retzlaff et al. 2008, in prep) offer us a superb op- 
portunity to constrain the volume density of very bright 
star- forming galaxies at z ~ 7—10. While these data 
are > 1 mag shallower than our NICMOS search data 
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TABLE 1 

Primary imaging data used for our z and J dropout searches. a 



5a Depth b 



Name 


Area 


2850 


J no 


#160 


K s 


ReF 


HDF-North Dickinson 


4.0 


27.6 


26.8 


26.8 


25.4 




1,2 




HDF-North Thompson 


0.8 


27.6 


27.8 


27.9 


25.4 




3,2 




HUDF Thompson 


5.8 


28.8 


27.4 


27.2 


25.8 




4,5 




HUDF Stiavelli 


0.7 


28.8 


27.9 


27.7 


25.8 




5.6 




HUDF-NICPAR1 


1.3 


28.4 


28.4 


28.2 






6.7 




HUDF-NICPAR2 


1.3 


28.4 


28.4 


28.2 






6,7 




GOODS Parallels 


9.3 


27.3 


26.9 


26.7 


~ 25 d 




[7] 


ISAAC v2.0 


136 


27.3 


~ 25. 5 d 


~ 25. d 


~ 25 d 


[8 


,9,10] 


MOIRCS GTO-2 


28 


27.3 


25.4 




25.4 


[2] 


MOIRCS GTO-1,3,4 


84 


27.3 


24.0 




24.2 




[2] 



a The layout of these search fields is illustrated in Figure 1. 

b 5<7 depths for ACS and NICMOS data given in terms of a 0.6"-diameter 
aperture and in terms of a ~ 1.0"-diameter aperture for the ground-based K s - 
band data. In contrast to the detection limits quoted in some of our previous 
work, here our detection limits have been corrected for the nominal light outside 
these apertures (assuming a point source). The detection limits without this 
correction are typically ~ 0.2 mag fainter. 

References: [1] Dickinson 1999, [2] Kajisawa et al. 2006, Ouchi et al. 2007, 
[3] Thompson et al. 1999, [4] Thompson et al. (2005), [5] Labbe et al. (2006), 
[6] Oesch et al. 2007, [7] Bouwens & Illingworth (2006), Riess et al. (2007), H. 
Teplitz et al. (2008, in prep) [8] B. Vandame et al. 2008, in prep; J. Retzlaff et 
al. 2007, in prep, [9] Mannucci et al. 2007, and [10] Stanway ct al. (2008). 

d The depth of the near-IR data over the CDF-South varies by ~ 0.2 - 0.4 
mag depending upon the observational conditions in which the ISAAC data 
were taken. 



considered above (§2.1), they cover > 5 times more area 
and thus are useful for constraining the volume density 
of the rarer, more luminous E/V-bright galaxies at z > 7. 
These data are most useful for searching for Zsscrband 
dropouts due to the greater depth of the ISAAC J-band 
data. While Mannucci et al. (2007) used an earlier reduc- 
tion of these data (vl.5) to search for bright z dropouts 
over the CDF-South GOODS field (see also Stanway et 
al. 2008), here we repeat this selection with our proce- 
dures. In doing so, we will be (1) able to more naturally 
include the search results over these fields with those ob- 
tained from our NICMOS fields and (2) be able to take 
advantage of the additional ISAAC data that have be- 
come available over the CDF-S since the completion of 
the Mannucci et al. (2007) work. The larger data set 
covers ~10% more area in the J and K s bands and 30% 
more area in the H band. We will not be using these im- 
ages to constrain the volume density of z ~ 9 J-dropouts 
because of the limited depth (~ 25 mag at her) of the 
ISAAC H and K s band data. 

We make use of the final (v2.0) reductions of the 
ISAAC data (B. Vandame et al. 2008, in prep; J. Retzlaff 
et al. 2008, in prep) for our searches. These reductions 
include ~ 143 arcmin 2 of deep J and K s band coverage 
and <~ 131 arcmin 2 of deep H band coverage. Of the area 
with deep J and K s coverage, only ~ 136 arcmin 2 has 
deep ACS GOODS coverage. This represents 4% more 
area than considered in the z-dropout search by Man- 
nucci et al. (2007) using the vl.5 ISAAC reductions over 
the CDF-South GOODS field. The FWHM of the PSF 
varies from <~ 0.35" to ~ 0.7", with a median value of 
~ 0.51". The depth of the J-band data is appropriately 
25.5 AB mag for a 4tr detection, while for the H and 
K s band data, the depth is ~ 24.9 AB mag. Of course, 
the actual depth varies by several tenths of a magnitude 



across the mosaic depending upon the precise exposure 
time of the ISAAC data and the FWHM of the PSF. As 
in the previous section, we used the same reductions of 
the ACS GOODS data as we used in our latest paper on 
the UV LF at z ~ 4, 5, and 6 (Bouwens et al. 2007). 

2.3. MOIRCS Search Fields 

We can strengthen our constraints on the prevalence 
of the bright star-forming galaxies at z > 7 by taking 
advantage of the very deep wide-area data taken over 
the HDF-North GOODS field. While there have been a 
number of independent attempts to obtain deep near-IR 
coverage over this field (e.g., the WIRC data over the 
HDF-North described by Erb et al. 2006), perhaps the 
best near-IR data set (i.e., deepest, widest area) is from 
the Subaru MOIRCS imaging observations (Kajisawa et 
al. 2006; Ouchi et al. 2007). This is due to the significant 
investment of telescope time (~ 17 nights) made by the 
MOIRCS GTO team in obtaining very deep J + K s im- 
ages during March, April, and May of 2006 (additional 
observations were made in March and April of 2007 but 
those are not yet public). There are also 2 nights of 
high quality K s band data available from a separate GO 
proposal. One of the main science goals for these obser- 
vations has been to search for z > 7 galaxies (e.g., Ouchi 
et al. 2007). 

The raw images were downloaded from the Subaru- 
Mitaka-Okayama-Kiso Archive (SMOKA: Baba et al. 
2002) and then divided into a number of discrete seg- 
ments, each segment containing <50 images taken one 
after another in a given filter. Images from each seg- 
ment were reduced and combined into a single image us- 
ing the iraf routines in MCSRED (I. Tanaka et al. 2008, 
in prep) . These reductions were then aligned against the 
ACS GOODS data and then stacked together to produce 
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two very deep near-IR image mosaics in the J and K a 
bands. These images were flux calibrated by matching 
sources with those in the KPNO images of the WFPC2 
HDF-North (Dickinson 1998) and in overlap regions be- 
tween adjacent MOIRCS exposures. Since the MOIRCS 
GTO team took their observations over 4 distinct point- 
ings, the final reductions cover ~ 112 arcmin 2 in total. 
The final reductions reach a 5a AB magnitude of 25.4 
in the J band and 25.4 AB mag in the K s band over 
the deepest MOIRCS pointing (GTO-2: - 28 arcmin 2 ), 
while the other three pointings (GTO-1,2,4) are some 
^1.3 mag shallower. The FWHM of the PSF in our re- 
ductions was <~ 0.5" in the K s band and <~ 0.6" in the 
J band. The final J-band image mosaic include some 
~56 hours of data, and the final K s band image mosaic 
includes some ~64 hours of Subaru data. 

3. SAMPLE CONSTRUCTION 

3.1. Catalog Construction 

Our procedures for generating PSF-matched cata- 
logues and selecting z and J-dropouts are already well- 
documented in several previous works (e.g., Bouwens et 
al. 2003; Bouwens et al. 2005; Bouwens & Illingworth 
2006; Bouwens et al. 2007). Nonetheless, we will briefly 
summarize the different steps in our processing here. All 
images were first PSF-matched to the NICMOS #160- 
band data (or J-band data in the case of our dropout 
searches over the ISAAC or MOIRCS data). We then 
ran SExtractor (Bertin & Arnouts 1996) in double-image 
mode to perform the object detection and photometry. 
We generated the detection image by taking the square 
root of x 2 image (Szalay et al. 1999) constructed from the 
Jno and -Hi6o band images for our z-dropout selection 
and from the i?i6o band image for our J-dropout selec- 
tion. Colors were measured in a small-scalable aperture 
using Kron-style (1980) photometry, where the Kron fac- 
tor was 1.2. These fluxes were then corrected up to total 
magnitudes using the light within a larger Kron (1980) 
aperture (adopting a Kron factor of 2.5). These latter 
corrections were made from the square root of the % 2 
image to improve the S/N. Figure 5 of Coe et al. (2006) 
provides a graphical description of a similar multi-stage 
procedure for measuring colors and total magnitudes. 
Aperture radii used for total magnitude measurements 
ranged from ~ 0.5—0.8 arcsec. An additional ~ 0.1 — 0.15 
mag correction was made to account for the light outside 
of these apertures and on the wings of the PSF (e.g., 
Thompson et al. 1998). 

3.2. Selection Criteria 

Our selection criterion for our z 850 -dropout search was 
((^850 - J\w)ab) > 0.8) A ((z 850 - Jiw)ab > 0.8 + 
0.4( Jno — -ffi6o) J 4s) as shown in Figure 2, where A repre- 
sents the logical AND operation. This criterion is iden- 
tical to the "less conservative" z-dropout selection crite- 
rion of Bouwens & Illingworth (2006). We opted to use 
this criterion rather the more stringent "conservative" 
criterion from that same work for two reasons. First, 
we realized that our "conservative" criterion may have 
been more stringent than necessary to reject contami- 
nants like T-dwarfs (low mass stars) from our samples. 
These contaminants can (in most cases) be rejected by 
a procedure we will describe below. Second, by adopt- 
ing a "less conservative" selection criteria, the effective 



selection volumes will be larger and therefore less sensi- 
tive to small inaccuracies in the assumed model galaxy 
population. 

For regions of the two GOODS fields with 
ISAAC/MOIRCS coverage, our z-dropout criterion 
was ((z 850 - J)ab > 1.0) A (J - K s < 1.2). This 
selection criterion is very similar to the criterion used 
over our NICMOS fields. For the typical star-forming 
galaxy at z ~ 7 (with a {/^-continuum slope (3 of —2) 
this criterion selects galaxies to a slightly lower-redshift 
limit z > 6.5 than our NICMOS z-dropout criterion 
(where z > 6.6). 

The criterion for our z ~ 9 Jno-dropout search was 
(Jno — H\6o)ab > 1-3. We elected to use a less re- 
strictive ( Jno — Hi 60 )ab color criterion than the (Jno — 
Hiqo)ab > 1-8 criterion used in our previous J-dropout 
search (Bouwens et al. 2005) because of the availabil- 
ity of deep optical data over almost all of our previous 
search fields to eliminate low redshift interlopers. Previ- 
ously ultra deep optical data were not available over the 
NICMOS parallels to the HUDF, and so we accordingly 
decided to select z dropouts using a more stringent J — H 
color cut. 

In computing the z — J and J — H colors that we used 
for our z and J dropout selections, we set the z- and 
J-band fluxes to their la upper limits in the case of non- 
detections. This was to minimize the extent to which 
faint sources in our search fields would scatter into our 
dropout samples. 

Sources in our z and J-dropout selections were also re- 
quired to be detected at the 4.5cr level in the #160 band 
within a 0.6" diameter aperture (to eliminate spurious 
sources) and undetected (< 2a) in the available -B435, 
Veoe, or *775 band imaging. Sources which were detected 
at the 1.5a level in more than one of the ACS £435, Vqoq, 
and Z775 bands were also rejected as low-redshift interlop- 
ers. Possible contamination by low mass stars (which are 
a concern for our z-dropout selection) are identified by 
looking for sources which showed faint detections (> 2a) 
in the z 850 -band in a 0.2"-diameter aperture and whose 
z 8 5o-band profile appeared to be consistent with a point 
source (i.e., the SExtractor stcllarity parameter > 0.6). 

All candidate dropouts were inspected in the avail- 
able IRAC and MIPS coverage over the two GOODS 
fields (Dickinson & GOODS team 2004) to ensure that 
none had particularly red 3.6/x — 5.8/i > 1 mag or 
3.6/i — 24/i > 2 mag colors. With the exception of 
the two sources which satisfy our J-dropout criterion 
(03:32:41.22, -27:44:00.8 and 03:32:32.45, -27:42:21.5), 
none of the sources which satisfied our other criteria were 
found to have such red infrared colors (or showed signif- 
icant detections at 5.8/x, 8.0/i, or 24/i). 

The two particularly red Jno dropouts (with i?i60 
magnitudes of 25.7 ± 0.2 and 24.5 ± 0.1, respectively) 
have #160 — [4.5/x] colors of ~ 2.5 mag and are even de- 
tected at 24/i in the available MIPS imaging over the 
GOODS fields (meaning that their apparent magnitudes 
at 24/im are < 21.5 mag). Since we do not expect galax- 
ies at z > 7 to be so bright at 24/n and otherwise so 
highly reddened (but see also discussion in Wiklind et 
al. 2007), we have excluded these sources as candidate 
z ~ 9 galaxies. 

3.3. Selection Results 
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V"110 AA 160/AB 

Fig. 2. — z — J and J — H color-color diagram illustrating our optical-IR 2-dropout selection window (orange region). The large black 
squares show the colors of z-dropouts in our selection (as well as one z-dropout found behind Abell 1689: Bradley et al. 2008) and their 
la errors (or lower limits). Sources which do not satisfy our z-dropout selection criteria are shown in cyan. The solid, dotted, and dashed 
blue lines show how the colors of starbursts with t/V-continuum slopes f3 of —2, —1, and 0, respectively, vary as a function of redshift. 
Typical f/V-continuum slopes /3 observed for star-forming galaxies at z ~ 5 — 6 are —2 (e.g., Stanway et al. 2005; Bouwens et al. 2006; Yan 
et al. 2005). The red lines show how these colors vary as a function of redshift for low-redshift galaxies (Coleman, Wu, & Weedman 1980). 
The green lines bracket the region in color-color space where we expect T dwarfs to lie (Knapp et al. 2004). The two red triangles indicate 
the colors of two bright T dwarfs found in our search (§3.3). While five of the z-dropout candidates in our selection nominally have z — J 
and J — H colors consistent with that of T-dwarfs, their z-band flux appears to be extended in all five cases (see also §3.2-§3.4), so these 
sources cannot be stars. 



Application of our selection criteria across all the fields 
given in Table 1 yielded 8 credible zgso-dropout can- 
didates, but no credible Jno-dropout candidates. All 
8 of our z-dropout candidates were found in our NIC- 
MOS search fields. The properties of each of the 8 Zs5o- 
dropout candidates are given in Table 2. Figure 3 pro- 
vides a postage stamp montage of the 8 candidates. We 
also included one candidate we identified in a separate 
search around massive lensing clusters (Bradley et al. 
2008; Bouwens et al. 2008). Five of the eight candidates 
were given in our two previously published Zgso-dropout 
searches (Bouwens et al. 2004c; Bouwens & Illingworth 
2006). Three of these candidates are new. Two of the 



three new candidates (UDF-3244-4727 and NICPAR1- 
3303-4111) were just below the flux thresholds used in 
our previous searches, but meet these thresholds now as 
a result of modest improvements in the S /N and overall 
quality of our NICMOS reductions. The third candi- 
date (NICPAR2-3308-5229) was found over the second 
NICMOS parallel to the HUDF and could be identified 
as a result of the much deeper optical (ACS) data now 
available over these fields to set strong constraints on the 
z — J colors. Our sample also includes one zsso-dropout 
(UDF-380-1125) which, while appearing in the our ini- 
tial z-dropout selection over the HUDF (Bouwens et al. 
2004c), did not appear in our more recent Bouwens & 
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Fig. 3. — V606i *775> z S50i JllOj an d ^160 band cutouts for 
all eight sources in our z ~ 7 zgso-dropout sample. The source 
presented in the bottom row is another very strong z ~ 7 z-dropout 
candidate found in an independent search around massive lensing 
clusters (Bradley et al. 2008; Bouwens et al. 2008). All the sources 
in our sample are detected at > 4.5er in both the Jno and -Higo 
bands, with the exception of UDF-3244-4727 and UDF-387-1125 
which is only detected at 3cr in the Jno band (we will be obtaining 
deeper Jno-band imaging data for these latter two sources as part 
of an approved HST program). All the sources in our sample show 
strong z — J breaks and arc completely undetected (< 2a) in the 
Veoe and 1775 bands (see §3.2). 

Illingworth (2006) selection because of the more conser- 
vative S/N threshold used in that work. 

There were two sources in our catalogues which while 
satisfying our z-dropout color criteria did not make 
it into our z-dropout selection because of their com- 
pactness in the z-band, which strongly suggested they 
were T dwarfs (03:32:25.10, -27:46:35.6 and 12:37:34.26, 
62:18:31.4). These sources had previously been iden- 
tified as possible T dwarfs in a number of other stud- 



ies (Bouwens & Illingworth 2006; Mannucci et al. 2007; 
Eylcs ct al. 2007). There were also 28 sources in our 
selection that while satisfying our z-dropout color cri- 
teria showed significant detections (> 2a) in the -B435, 
V6O61 ° r «775 bands and thus were excluded from our 
z-dropout selection. One possible z-dropout candidate 
over the near-IR ISAAC data was found at 03:32:24.74, 
—27:55:15.8. This source seemed to be a plausible z ~ 7 
candidate, with an apparent 5a detection in the J band, 
an apparent J-band magnitude of 24.8±0.2 AB mag, and 
no detection (< la) in any of the optical imaging data. 
However, since this source is undetected (Hab > 26.0 [3a 
limit] in a l"-diameter aperture) in even deeper NICMOS 
i?i6o-band data over this field, it seems highly likely that 
this source is spurious (or perhaps a SNe: see §3.4). 

In Appendix D, we discuss several sources which ap- 
peared in previous z > 7 dropout samples of ours but do 
not appear in the current selections. 

3.4. Contamination 

We give particular attention to contamination, since 
it can be a significant issue in searches for very faint, 
high-redshift galaxies. This discussion is only relevant 
for our Zsso-dropout selection since our Jno-dropout se- 
lection contains no credible candidates. There are four 
types of contaminants that are of particular concern: (i) 
T-dwarfs, (ii) supernovae, (Hi) lower-redshift galaxies 
which may move into our Zsso-dropout selection through 
photometric scatter, and (iv) spurious sources. 

Possible Contamination from T dwarfs: T dwarfs are 
low-mass brown dwarfs, which while most well known 
and studied in the solar neighborhood (e.g., Knapp et 
al. 2004; Burgasser et al. 2006), are sometimes observed 
at great distances and therefore can be quite faint (i.e., 
Jab ~ 26: Ryan et al. 2005; Bouwens & Illingworth 
2006). Since T dwarfs occupy a position in color-color 
space very similar to that of high-redshift galaxies (Fig- 
ure 2), they can therefore be difficult to distinguish from 
star-forming galaxies at z ~ 7 on the basis their photom- 
etry alone. This leaves us with the size or morphology of 
these sources as essentially our only means of distinguish- 
ing these sources from z <~ 7 galaxies. This enterprise is 
made more difficult because the T dwarfs are brightest 
in exactly the same data (near-IR) where the resolution 
is the poorest. Fortunately, the available ACS data have 
sufficient depth and resolution that any T dwarfs in our 
data should show up as > 5a point sources. This fol- 
lows from the fact that all candidate dropouts in our 
sample are detected at > 4.5er, T dwarfs have z ~ J col- 
ors of ~ 1 — 2 mag, and the ACS zsso-band data are at 
least ~ 2 mag deeper in a 0.2"-diameter aperture (ACS 
resolution) than our near-IR data is in a 0.6"-diameter 
aperture (near-IR resolution: see Table l). 1 Since none 
of the z-dropouts in our selection show > 4cr-pointlikc 
detections in the available zsso-band coverage, this sug- 
gests that these sources are not T dwarfs. 

Possible Contamination from Supernovae: Since the 
NICMOS JH data were largely taken at quite different 
times than the ACS BViz data, it is quite possible that 
highly time-variable sources like supernovae could show 

1 The 5<t depths of the ACS data in a 0.2"-diameter aperture is 
~ 1.2 mag deeper than that given in Table 1 (for a 0.6"-diamctcr 
aperture). 
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TABLE 2 

Z850-DROPOUT CANDIDATES. 



Object ID 


R.A. 


Dec 


#160 


2850 


— Jno 


</no 




#160 


#160 — K s 


ReP 


HDFN-3654-1216 


12:36:54.12 


62:12:16.2 


26.0 


±0.1 


1.1 


±0.3 


0.4 


± 


0.3 


0.0 


±0.3 


[2] 


UDF-640-1417 


03:32:42.56 


-27:46:56.6 


26.2 


±0.1 


1.3 


±0.3 


0.5 


± 


0.2 


0.5 


±0.3 


[1,2,3,4,5] 


NICPAR2-3308-5229 b 


03:33:08.29 


-27:52:29.2 


26.7 


±0.1 


> 


2.2 C 


0.6 


± 


0.2 








CDFS-3225-4627 


03:32:25.22 


-27:46:26.7 


26.7 


± 0.2 


1.4 


±0.5 


0.2 


± 


0.3 


< 


0.2 C 


[2] 


UDF-983-964 


03:32:38.80 


-27:47:07.2 


26.9 


±0.2 


> 


3.2 C 


0.1 


± 


0.3 


0.0 


±0.8 


[1,2,4,5] 


UDF-387-1125 


03:32:42.56 


-27:47:31.4 


27.1 


±0.2 


1.3 


±0.5 


0.8 


± 


0.4 


< 


0.0 C 


[1,5] 


UDF-3244-4727 


03:32:44.02 


-27:47:27.3 


27.3 


±0.2 


> 


2.6 C 


0.5 


± 


0.4 


< 


0.1 c 




NICPARl-3303-4111 d 


03:33:03.81 


-27:41:12.1 


27.8 


±0.1 


> 


1.5° 


0.4 


± 


0.2 








A1689-zDl e 


13:11:29.73 


-01:19:20.9 


24.7 ±0.1 


> 2.2 C 


0.6 ±0.2 






[6] 



a Rcfcrcnccs: [1] Bouwcns ct al. (2004c), [2] Bouwens & Illingworth (2006), [3] Yan & Windhorst (2004), [4] Coc ct al. 
(2006), [5] Labbc ct al. (2006), [6] Bradley et al. (2008) 

b There appears to be some flux in the Veoe band very close to the position of this source (~ 0.1"). If this flux is 
associated with this source and not a chance superposition, it would rule out this source being a z > 7 galaxy. The 
existence of a small amount of flux in the Veoe band would be somewhat unexpected given that there is no detectable 
flux in the 1775 and zgso bands and the fact that these exposures reach similar flux limits. On balance, we consider this 
to be a likely z ~ 7 candidate. 

c Upper and lower limits on the measured colors are the lcr limits. 

d Bccause this source (detected at ~ 5 — 6<r) was found over the NICMOS parallels to the HUDF and the noise properties 
of these data may exhibit significant non-Gaussian characteristics due to the small number of unique dither positions 
(see §2.1), we took special care in assessing the reality of this source. We considered all the available NICMOS data at 
4 highly distinct dither positions (separated by ~ 2 — 3" each) and in both the J110 and Hi6o-band data separately. 
The candidate was evident at > 2 — 3<r significance in 7 of the 8 of the stacks, consistent with what we might expect 
if the candidate were a real source. This source is also present in an independent reduction of these data (Ocsch et al. 
2008). 

°This source was identified in a separate search for z ~ 7 z-dropouts behind massive galaxy clusters (Bradley ct al. 
2008; Bouwens et al. 2008) and is not used here for deriving the UV LF at z ~ 7. This source is substantially magnified 
by gravitational lensing from Abell 1689 (/1 ~ 9.3). Its delensed /fi6o-band magnitude is ~ 27.1 AB mag (Bradley et 
al. 2008). 



up as dropouts in our selections. To estimate the approx- 
imate frequency at which we might expect such sources 
to show up in our searches, we adopted the GOODS su- 
pernova searches (Riess et al. 2004; Strolger et al. 2004) 
MS M baseline. In their searches, they identified 42 su- 
pernovae in 8 independent searches over one of the two 
~ 160 arcmin 2 GOODS fields (15 ACS WFC pointings), 
or ~ 1/3 supernovae per ~ 11 arcmin 2 ACS WFC field 
per ~ 40 day period. Of these supernovae, only ~ 40% 
of the high-redshift supernovae identified by Riess et al. 
(2004) and Strolger et al. (2004) appear to be well enough 
separated from their host galaxies to be identified as a 
separate source. Having calculated the approximate rate 
at which supernovae appear at random in deep fields, all 
that remains is to calculate the time period over which 
we might expect these supernovae to show up as possi- 
ble sources in our search. While theoretically we might 
expect high-redshift SNe to be bright enough to show up 
in our survey data for ~ 100 days (if we take the typ- 
ical light curves found in e.g., Strolger et al. 2004), all 
z-dropouts in our sample are within ~ 1 mag of the flux 
limit of each survey field. Typical SNe (at z ~ 1) would 
not remain at these magnitudes for more than ~ 40 
days, so the incidence of SNe in our probe should be 
- (l/3)(l/llarcmin 2 )(0.40) ~ 0.012 arcmin" 2 for data 
taken at different epochs. Multiplying this by the ^16 
arcmin 2 of NICMOS data in our search taken at signif- 
icantly different times than the ACS data (the HUDF 
NICMOS data were taken at a similar time to the ACS 
HUDF data), we would expect that only ^0.2 supernovae 
should contaminate our high redshift selections. While 
this number is just 3% of the 8 sources in our z-dropout 
selection and therefore likely of minimal importance for 



the present LF determinations, it does indicate that con- 
tamination of high-redshift selections by supernovae is a 
real possibility. This may be a particularly important 
concern in searches for bright z > 7 galaxies due to the 
inherent rarity of those sources. 

Possible Contamination Through Photometric Scatter: 
One possible source of contamination are low redshift 
sources scattering into our color selection simply because 
of noise. To assess the importance of this effect, we ran 
a number of Monte-Carlo simulations where we started 
with a low-redshift galaxy population and then added 
noise to the fluxes in the individual bands to see how 
many of the sources would satisfy our z 850 -dropout se- 
lection. Clearly, the most critical aspect of these simula- 
tions is ensuring that the properties of the input galaxy 
population (to which noise is added) are accurate. We 
set up the input population as follows: (1) The input 
Hi60,AB magnitudes and photometric errors are mod- 
elled using the actual distribution of magnitudes and er- 
rors found in the data (to make our estimates as model 
independent as possible and because essentially all faint 
sources in our fields will be at z < 6.5) and (2) The colors 
(i.e., optical-i/ 160 and J 110 — i?i6o) of the input popula- 
tion arc assumed to be distributed in the same way as is 
observed for galaxies in the range 24.5 < -ffi6o ab < 26 in 
the HUDF (ACS+NICMOS coverage). We used galax- 
ies in this magnitude range because of their higher S/N 
colors and because the multivariate color distribution is 
similar (albeit slightly redder) to galaxies at 7Ji6o-band 
magnitudes > 26 AB mag. This magnitude range also 
does not include any sources which are candidate z ~ 7 
galaxies (at least in the HUDF). 

Based upon the described input population, we ran 
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10 separate photometric scattering experiments for each 
search held. We estimated that there would be 0.0, 
0.3, 0.2, and 0.3 low-redshift contaminants, respectively, 
for our z-dropout selections over the HUDF, NICPAR, 
HDFN, and GOODS parallel fields. In total, - 0.8 con- 
taminants are expected from all our fields. We note that 
the approach taken here to estimate contamination from 
photometric scatter is analogous to the one described in 
Bouwens et al. (2004c) or Appendix D.4.2 of Bouwens et 
al. (2006). 

Possible Contamination from Spurious Sources: Fi- 
nally, to estimate the number of contaminants from spu- 
rious sources, we repeated our selection procedure on the 
negative images (e.g., Dickinson et al. 2004; Yan & Wind- 
horst 2004; Bouwens et al. 2007) and found no sources 
that satisfied our selection. This suggests that the spu- 
rious fraction is negligible. We note that we would not 
expect contamination from spurious sources to be a prob- 
lem since essentially all of the sources in our sample are 
detected at > 4.5ct in both the Jno and i?i60 bands. The 
only notable exceptions to this are for the two sources in 
our sample, UDF-387-1125 and UDF-3244-4727, where 
Jno-band detections are still only at the 3<r level. How- 
ever, in both cases, we are attempting to improve the 
detection significance (as well as the S/N of the z — J 
and J — H colors) by obtaining an additional 9-12 orbits 
of Jno-band data on both sources (likely improving their 
detection significance to the 4.5cr level). 

In summary, the overall level of contamination should 
be fairly low in general. We would expect ~ 0.8 galaxy 
to enter our sample from photometric scatter, <~ 0.2 SNc 
contaminants to enter our sample, and there to be no 
contamination from T dwarfs or spurious sources. Over- 
all, this suggests an overall contamination level of ~ 12%. 

4. DETERMINATIONS OF THE UVLFATZ~7-10 

We now utilize our search results to determine the rest- 
frame UV LF at z ~ 7. We will consider both LFs using 
both stepwise and Schechter parametrizations. The ad- 
vantage of the stepwise determinations is that they give 
a more model-independent measure of the constraints we 
have on the LF, while the Schechter determinations are 
more amenable to direct interpretation and comparison 
with other determinations in the literature. 

As in other methodologies for determining the LF, we 
maximize the probability with which our model LFs are 
able to reproduce our search results. Since we do not 
know the precise redshifts for our sources, our goal will 
simply be to reproduce the surface density of z-dropout 
candidates in each of our search fields as a function of 
apparent magnitude. This is similar to the modelling 
we performed in our previous paper on the UV LFs at 
z ~ 4, z ~ 5, and z ~ 6 (Bouwens et al. 2007), where we 
attempted to reproduce the surface density of B, V, and 
i dropouts in deep ACS helds. 

To maximize the likelihood with which our model LFs 
reproduce our search results, we need to be able to com- 
pute the expected surface density of dropouts as a func- 
tion of magnitude. We use a similar formalism to the one 
we used in our previous paper on the LFs at z ~ 4 — 6: 

^k4>kV m ,k = N m (1) 

where N m is the surface density of galaxies in some search 
field with magnitude m, 4>u is the volume density of 



galaxies with absolute magnitude k, and V m .k is the ef- 
fective selection volume for which galaxies with absolute 
magnitude k will both satisfy our dropout selection cri- 
teria and be observed to have an apparent magnitude m. 
The width of the magnitude bins in N m , fa, and V m ^ 
will depend upon the LF parameterization we consider 
(stepwise vs. Schechter). 

We computed the effective volumes V m ,k for our z 
and J dropout selections by artificially redshifting B- 
dropouts from the HUDF across the selection windows 
of our higher rcdshift samples, adding the sources to the 
real data, and then repeating our selection procedure. 
We used our well-tested cloning software (Bouwens et 
al. 1998a, b; Bouwens et al. 2003) to handle the artificial 
redshifting of sources from z~4toz~6 - 10. The 
size of the sources were scaled as (1 + z) -11 (indepen- 
dent of luminosity) to match the size trends observed for 
star-forming galaxies at z > 2 (Bouwens et al. 2006; see 
also Ferguson et al. 2004 and Bouwens et al. 2004c). The 
rcdshifted ( "cloned" ) sources were smoothed so that their 
effective PSFs matched that seen in the data. The mean 
[/V-continuum slope (3 was set to —2 to agree with that 
observed in L* =3 star-forming galaxies at z ~ 6 (Stan- 
way et al. 2005; Bouwens et al. 2006; Yan et al. 2005) 
while the la dispersion in (3 was set to 0.5, which is sim- 
ilar to what is found at z ~ 4 — 6 (e.g., Bouwens et al. 
2007; R.J. Bouwens et al. 2008, in prep). Instead of us- 
ing simple power laws to represent model SEDs of given 
UV continuum slope [3, we elected to use 10 8 -yr con- 
tinuous star-formation models (Bruzual & Chariot 2003) 
where the dust extinction (Calzetti et al. 1994) is var- 
ied to reproduce the model slopes. This should provide 
for a slightly more realistic representation of the SEDs 
of star-forming galaxies at z > 7 than can be obtained 
from simple power law spectra. 

To evaluate the likelihood with which the surface den- 
sities computed above would reproduce that found in 
the observations, we assumed simple Poissonian statis- 
tics and compared the observed and expected numbers 
bin by bin. In detail, our expression for the likelihood C 
is 

£ = n«e-^%^ (2) 

where N obs _ij is the observed number of sources in search 
field i and magnitude interval j, where N exPi ij is the 
expected number of sources in search field i and magni- 
tude interval j, and where ILij is the product symbol. 
We elected to use the above procedure instead of a more 
conventional procedure like STY79 (Sandage, Tammann, 
& Yahil 1979) or SWML (Efstathiou et al. 1988: which 
determine the shape of the LF independent of the nor- 
malization) to take advantage of our search constraints 
on the volume density of bright z > 7 galaxies from our 
wide-area ground-based data. Maximum likelihood pro- 
cedures (like SWML or STY79) do not take advantage 
of these null search results - simply because there are 
no sources to use for estimating the shape of the LF. 
While the maximum likelihood fit results from Eq. 2 will 
be affected by large-scale structure variations across our 
search fields, we will compute the uncertainties which re- 
sult from these variations in §4.1, §4.2, and Appendix A 
and include them in our overall error budget. 
In determining the rest-frame UV LF at z ~ 7 and 
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z ~ 9 from our z and J dropout selections, we will adopt 
z <~ 7.3 and z ~ 9.0, respectively, as the fiducial redshifts 
for these selections and the derived LFs. In §4.2, we will 
justify the use of these mean redshifts using the effective 
volume simulations just described (see also Bouwens et 
al. 2004; Bouwens et al. 2005; Bouwens & Illingworth 
2006). 

4.1. Stepwise Determinations 

We begin with a determination of the LF at z ~ 7 
and z ~ 9 in stepwise form. The advantage of step- 
wise parametrizations is that they allow us to determine 
the overall prevalence of galaxies in a way that is very 
model independent and just a simple function of lumi- 
nosity. The width of the bins in the LF are 0.8 mag and 
was chosen as a compromise between the need to have a 
few z ~ 7 galaxies per bin and the desire to retain a rea- 
sonable resolution in luminosity for our derived LFs. We 
use a similar bin size for the stepwise upper limits on our 
z ~ 9 J-dropout LF in order to make these limits more 
easily interpretable. The stepwise LF fa is then derived 
using Eq. 1. The results are presented in Table 3. 

Perhaps the largest uncertainties in these LF determi- 
nations arise from small number statistics (Poissonian er- 
rors). However, in addition to the simple number statis- 
tics, there are also sizeable uncertainties that come from 
large-scale structure variations as well as our effective 
volume estimates. Let us first consider the uncertainties 
coming from large-scale structure variations. Adopting 
a pencil beam geometry, redshift selection window with 
width Az = 1.5, search area of 100 arcmin 2 , and bias of 
7 (which is appropriate for sources with volume densities 
of ~ 10" 3 - 5 Mpc -3 : e.g., Mo & White 1996; Somervillc 
et al. 2004; Trenti & Stiavelli 2008), we estimate a lcr 
RMS uncertainty of ^30% in the </>fc's due to field-to-field 
variations (see also estimates in Bouwens & Illingworth 
2006). At z ~ 9, these uncertainties are ~ 20% assum- 
ing a redshift selection window with width Az = 2 (see 
Bouwens et al. 2005). 

Uncertainties in our effective volume estimates derive 
primarily from our imperfect knowledge of the size (or 
surface brightness) distribution of star-forming galaxies 
at z > 7. 2 Fortunately, the mean size of star-forming 
galaxies at z > 4 show a good correlation with redshift 
(i.e., mean half-light radius cx (l + z) _11±0,3 for fixed lu- 
minosity: Bouwens et al. 2006) and we can make a rea- 
sonable estimate for what the size (surface brightness) 
distribution of galaxies is at z > 7. Nonetheless, this 
distribution is at least as uncertain as the error on the 
size-redshift scaling. Propagating the error on this scal- 
ing into the size distributions assumed in our effective 
volume estimates, we estimate an RMS uncertainty of 
17% in the selection volume at z ~ 7 and 15% at z ~ 9 
due to the uncertainties in the size (surface brightness) 
distribution. Together the size and large-scale structure 
uncertainties add an uncertainty of 34% and 25% RMS 
to each bin of the rest-frame UV LF at z ~ 7 and z ~ 9, 
respectively. These uncertainties have been added in 
quadrature with those deriving from the small number 
statistics. They are given in Table 3. 

2 Uncertainties in the UV color distribution also contribute to 
the overall error budget for our effective volume estimates, but they 
are smaller in general (e.g., see §5.2). 
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TABLE 3 

Stepwise Constraints on the 
rest-frame uv lf at 2 ~ 7 and 
2-9 (§4.1). 

M uv>AB b <t> k (Mpc- 3 mag" 1 ) 

z-dropouts (z ~ 7) 
-21.55 < 0.000003 a 

-20.84 0.00005lg;K 
-20.04 
-19.24 

J-dropouts (z ~ 9) 
-21.95 < 0.000019 a 

-21.15 < 0.000023 a 

-20.35 < 0.000082 a 

-19.55 < 0.000687 a 



a Uppcr limits here are lcr (68% 
confidence) . 

b The effective rest-frame wave- 
length is ~ 1900A for our 2-dropout 
selection and ~ 1500A for our J- 
dropout selection. 

These LFs are also presented in Figure 4 with the ma- 
genta points for our z ~ 7 LF and black downward arrows 
for the constraints on the z ~ 9 LF. A comparison with 
previous determinations at z ~ 4, z ~ 5, and z <~ 6 from 
Bouwens et al. (2007) is also included on this figure for 
context. Though the error bars for individual points in 
the LF at z ^ 7 are still quite sizeable, there is strong ev- 
idence that the UV LF at z ~ 7 is different from the UV 
LF at z ~ 6 (99% confidence) and thus there is evolution 
from z <~ 7 to z ~ 6. We determined this confidence 
level by finding the value of M* and <f>* which minimizes 
the total \ 2 evaluated for our i and z dropout LFs and 
then looking at the probability of obtaining the resul- 
tant reduced-^ 2 purely by chance. This conclusion was 
already drawn by Bouwens & Illingworth (2006) on the 
basis of a smaller but very similar selection of galaxies. 

4.2. Schechter Determinations 

We will now attempt to express the results of our 
search for z > 7 galaxies using a Schechter parameteriza- 
tion. The Schechter parametrization is convenient since 
it is much more amenable to interpretation than stepwise 
LFs are. Of course, it is not at all clear from the step- 
wise LFs derived in §4.1 (particularly given the sizeable 
observational uncertainties) that the UV LF at z > 7 is 
well described by a Schechter function (see discussion in 
§5.5). 

As with our stepwise determinations, we will calculate 
the expected surface density of dropouts given a model 
LF by using Eq. 1 and expressing the Schechter function 
in stepwise form. For convenience, we have decided to 
bin the surface density of galaxies in magnitude intervals 
of width 0.1 mag. Use of substantially finer bins does not 
have a noticeable effect on the results. Because of the size 
of current z ~ 7—10 samples and limited luminosities 
(< —19 AB mag) to which we can probe, we cannot hope 
to obtain very strong constraints on the faint-end slope 
of the LF at z > 7 and therefore it makes sense for us to 
fix it to some fiducial value. We will adopt —1.74, which 
is the faint-end slope of the UV LF at z ~ 6 determined 
by Bouwens et al. (2007) using the HUDF and a large 



11 




'"UV.AB 

Fig. 4. — Determinations of the rest-frame UV luminosity func- 
tion (LF) at z ~ 7 using both a Schechter parameterization (ma- 
genta line) and in stepwise form (magenta circles with lcr error 
bars). Note that the stepwise and Schechter determinations of the 
LF are determined separately (i.e., our Schechter LF fits are not 
obtained through fits to our stepwise LFs). The lines are not fits 
to the points. The lcr upper limits on the bright end of the UV 
LF at z ~ 7 and at z ~ 9 are shown with the downward arrows 
in magenta and black, respectively. For context, we have included 
the rest-frame UV LFs determined by B07 at z ~ 4 (red symbols), 
z ~ 5 (green symbols), and z ~ 6 (red symbols). 



TABLE 4 

Comparison of the Schechter Parameters for the 
rest-frame uv lfs at z ~ 7 and 2 ~ 9 with determinations 
at 2 ~ 4, 2 ~ 5, and z ~ 6 from bouwens et al. (2007). 



Dropout 
Sample 


< z > 




<j>* (10 -3 
Mpc -3 ) 


a 


z 
J 


7.3 
9.0 


-19.8 ±0.4 
> -19.6 C 


(l.l)b 


(-1.74) b 
(-1.74) b 


B 
V 

i 


3.8 
5.0 
5.9 


-20.98 ± 0.10 
-20.64 ±0.13 
-20.24 ±0.19 


1.3 ±0.2 
1.0 ±0.3 

1 4+0-6 


-1.73 ±0.05 
-1.66 ±0.09 
-1.74 ±0.16 



a Values of M* v are at 1600 A for the Bouwens et al. (2007) B 
and V-dropout LFs, at ~ 1350 A for the Bouwens et al. (2007) i- 
dropout LF, at ~ 1900 A for our 2-dropout LF, and at ~ 1500 A 
for our constraints on the J-dropout LF. Since z ~ 6 galaxies are 
blue (j3 ~ —2: Stanway et al. 2005; B06), we expect the value of 
My V at z ~ 6 to be very similar (< 0.1 mag) at 1600 A to its value 
at 1350 A. Similarly, we expect the value of M^ v at z ~ 7 and 
z ~ 9 to be fairly similar at ~ 1600A to the values determined at 
~ 1900A and ~ 1500A, respectively. 

b Fixed in our fits (§4.2) 

Tower limits here are lcr (68% confidence). 



number of deep ACS fields. Later we will investigate the 
sensitivity of our results to the value we assume for the 
faint-end slope. 

Our best-fit Schechter parameters at z ~ 7 for our 
z-dropout selection are M* AB = —19.8 ± 0.4 mag and 
<p* = 0.00111^07 Mpc- 3 for a fixed faint-end slope 
a = —1.74. The 68% and 95% likelihood contours for 
these parameters are given in Figure 5 and compared 
with our previous determinations from our B, V, and 
i-dropout selections at z ~ 4, z ~ 5, and z ~ 6, respec- 




FlG. 5. — 68% and 95% likelihood intervals on the Schechter 
parameters M* and <f>* for the rest-frame UV luminosity function 
(LF) at z ~ 7 (magenta contours) and z ~ 9 (black contours) de- 
rived from our z and J-dropout search results, respectively (where 
the left and right black lines represent the 95% and 68% confidence 
intervals, respectively, on the z ~ 9 LF). Also shown are the like- 
lihood contours Bouwens et al. (2007) derived for the UV LF at 
z ~ 4 (blue contours), z ~ 5 (green contours), and z ~ 6 (red 
contours). These contours are expressed in terms of the Schechter 
parameters c6* and M*. A faint-end slope a of —1.74 is assumed 
for the determinations at z ~ 7 and 2 ~ 9 (to match that seen at 
z ~ 6: Bouwens et al. 2007), but the exact value assumed for a 
makes little difference to the derived contours for <f>* and M* (see 
§4.2). Not surprisingly, our constraints on the Schechter parame- 
ters at 2 ~ 7 are considerably worse than at z ~ 4 — 6, reflecting 
the substantially smaller size of samples at z > 7. Nonetheless, it 
seems clear that the UV LF shows significant and consistent evo- 
lution towards brighter values of M* with cosmic time (by ~ 1.2 
mag from 2 ~ 7 to 2 ~ 4). There is little evidence for evolution in 
the volume density c6* . Note that our constraints on the z ~ 9 LF 
rule out no-evolution from 2~9to2~7at~ 80% confidence. 
No evolution from z ~ 7 to z ~ 6 is ruled out at 99% confidence. 

tively (Bouwens et al. 2007). The best-fit values are also 
given in Table 4. Large-scale structure uncertainties re- 
sulting from ficld-to-ficld variations were estimated using 
Monte-Carlo simulations (Appendix A) and incorporated 
into the uncertainties quoted above. While the best-fit 
value for cf>* is very similar to that found at z ~ 4 for the 
Bouwens et al. (2007) 73-dropout selections, the best-fit 
value for Mj} v is 1.2 ± 0.4 mag fainter than the value of 
M(j V (=-20.98 ± 0.07 mag) found at z ~ 4 by Bouwens 
et al. (2007) and 0.4 ± 0.4 mag fainter than the value of 
M* v (= -20.24 ±0.19 mag) found at z ~ 6 by Bouwens 
et al. (2007). This suggests that the brightening we ob- 
serve in My V from z ~ 6 to z ~ 4 (Bouwens et al. 2006; 
Bouwens et al. 2007) is also seen from z ~ 7. Of course, 
we must admit that we are somewhat surprised that our 
best-fit Schechter parameters are in such excellent agree- 
ment with an extrapolation of lower redshift trends! It 
would suggest that our z-dropout sample may largely be 
made up of star-forming galaxies at z ~ 7 as we have 
argued in §3.3 and §3.4 (i.e., the number of low-rcdshift 
interlopers is small) and that the effective volumes we 
have estimated for this sample are reasonably accurate 
(see also discussion in Appendix B). 
Given the small size of current z-dropout samples, it 
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Fig. 6. — Number of z-dropouts in our search fields (histogram) 
versus that expected from our best-fit LF (red lines). Panels cor- 
respond to searches over the two NICMOS parallels to the HUDF 
(Bouwens et al. 2005; Oesch et al. 2007), the NICMOS HUDF 
Thompson & Stiavelli fields (Thompson et al. 2005; Oesch et 
al. 2007), the HDF-North Dickinson (1999) field, the HDF-North 
Thompson et al. (1999) field, the CDF-South GOODS ISAAC 
field (B. Vandame et al. 2008, in prep; J. Retzlaff et al. 2008, in 
prep: §2.2), the HDF-North GOODS MOIRCS fields (Kajisawa et 
al. 2006: §2.3) and various NICMOS parallels over the GOODS 
fields (see Table 1 and Figure 1). Binning is done on 0.5 mag in- 
tervals. The green line shows the expected counts if we force our 
best-fit UV LF at z ~ 7 to have the same value of M^j V and a 
as we determined at z ~ 4 and only allow <fi* to vary to best fit 
the data. Clearly, the best-fit LF obtained in this manner (green 
lines) predicts too many galaxies at bright magnitudes than are ob- 
served (e.g., ~ 2 — 3 galaxies predicted at Jab, Hl6Q,AB ~ 25 — 26 
vs. the observed) and far too few galaxies at fainter magnitudes 
(-^160, AB ~ 27 — 28: particularly in the HUDF panel where the 
observed numbers are > 4x larger). This demonstrates that even 
in the very modest z-dropout selections which are available there is 
reasonable (~ 3<r) evidence that the typical luminosity of galaxies 
at z ~ 7 is much fainter than it is at z ~ 3 — 4. Of course, despite 
our ability to draw a few simple conclusions like this, there are still 
significant uncertainties in deriving the UV LF from the observed 
counts. 



may seem surprising that we are able to obtain any con- 
straint at all on the shape on the UV LF at z ~ 7. For- 
tunately, the large luminosity range over which we have 
constraints on the surface density of dropouts (i.e., from 
25 AB mag to 28 AB mag) largely makes up for what we 
lack in statistics. These constraints can be helpful, even 
brightward of 26.0 AB mag where our z-dropout sample 
contains no sources, because of the large search areas we 
consider (^248 arcmin 2 ). To illustrate how our search 
results help us constrain the overall shape of the UV LF, 
we show in Figure 6 the predicted counts for each of our 



TABLE 5 

Best-fit schechter parameters for the UV LFs at 
Z ~ 7 AND z ~ 9 WITH DIFFERENT assumed faint-end 
SLOPES. 



Dropout 






<t>* (10" 3 






Sample 


< z > 


lu uv 


Mpc~ 3 ) 




a 


z 


7.3 


-19.8 ±0.4 


J - 1 -0.7 


( 


-1.74) b 


z 


7.3 


-19.6 ±0.4 


i 7+2.6 
s -'—\.l 




-1.4) b 


z 


7.3 


-19.9 ±0.4 






-2.0) b 


J 


9.0 


> -19. 6 C 


(l.l) b 


( 


-1.74) b 


J 


9.0 


> -19. 5 C 


(l.l) b 




-1.4) b 


J 


9.0 


> -19. 7 C 


(l.l) b 




-2.0) b 


J 


9.0 


(-21.0) b 


< 0.06 c 


( 


-1.74) b 


J 


9.0 


(-21.0) b 


< 0.05 c 




-1.4) b 


J 


9.0 


(-21.0) b 


< 0.07 c 




-2.0) b 



a See remarks in Table 4. 
b Fixed in our fits (§4.2) 

Tower limits here are la (68% confidence). 



search fields from our best-fit UV LF (i.e., Mfj V — —19.8 
mag, 4> = 0.0011 Mpc -3 ) in red and compare it against 
the best-fit model where My V shows no evolution from 



4 (green lines: M ( 



uv 



-21, 



0.00012 Mpc~ 



The observed surface densities of dropouts are shown 
with the solid histogram. As is apparent from this fig- 
ure, the model with the bright value of My V predicts 
too many bright sources (e.g., ISAAC/MOIRCS panel) 
and too few faint sources (e.g., HUDF panel), whereas 
the model with a fainter value of M^ v provides a better 
match to the observations. The significance of this result, 
however, is clearly only modest (3<x) at present, given the 
limited statistics. Better statistics would also be helpful 
for verifying that our treatment of the selection effects in 
each search field is reasonably accurate. 

For our z ~ 9 J-dropout selection, we cannot derive 
a set of best-fit Schechter parameters due to the lack of 
any dropout candidates in our sample. Nonetheless, it 
does make sense for us to set constraints on Mfj V and <f>* 
assuming no evolution in the other Schechter parameters 
from z ~ 4. For example, supposing there is no evolution 
in My V (i.e., M'(j V = —21), we can set a la upper limit 
of < 0.00006 Mpc -3 (~ 20 x lower than at z ~ 4) on 
<p* . Similarly, assuming there is no evolution in <fi* (i.e., 
4>* = 0.0011 Mpc~ 3 ), we can set a la lower limit of —19.6 
AB mag (1.4 mag fainter than at z ~ 3) on M^ g . The 
joint constraints on the value of <jf and M^j v are shown 
in Figure 5 with the black lines indicating the 68% and 
95% confidence limits. 

In constraining the Schechter parameters at z ~ 7 and 
z ~ 9 from our z and J-dropout searches, we assumed 
that the faint-end slope a was equal to the same value 
we had previously determined at z ~ 6, i.e., a = —1.74. 
To see how these constraints might depend upon the as- 
sumed faint-end slope a, we repeated the determinations 
assuming a faint-end slope a of —1.4 and —2.0. The re- 
sults are presented in Table 5, and it is clear that the 
derived values for Myy and (jf only show a mild de- 
pendence on the assumed faint-end slope a, particularly 
when compared to the uncertainties on these two param- 
eters from other error sources (e.g., from small number 
statistics). We can therefore essentially ignore the effect 
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Fig. 7. — The redshift distributions predicted for our z and 
J-dropout selections (shown as the magenta and black lines, re- 
spectively) using our derived LFs and the simulation results pre- 
sented at the beginning of §4. The z-dropout predictions assume 
Muv,AB = -19.8 and a = -1.74 (Table 4) while the J-dropout 
prediction assume Mjjv,ab = —19.1 and a = —1.74 (derived from 
the extrapolations derived in §5.3). The mean redshifts expected 
for our z-dropout and J-dropout selections are 7.3 and 9.0, re- 
spectively. These mean redshifts are somewhat lower than that of 
previous selections we had considered (e.g., Bouwens et al. 2004c; 
Bouwens et al. 2005; Bouwens & Illingworth 2006) because of our 
use of slightly more inclusive selection criteria in this work. 

of the uncertain faint-end slope on our quoted constraints 
for both M^y and <fi* . 

To put the LF fit results and dropout selections in 
context, it seems helpful to use our best-fit LFs to cal- 
culate an expected redshift distribution for our z and 
J dropout selections. We can take advantage of the 
same simulations described at the beginning of §4 to 
calculate these redshift distributions. For concreteness, 
MijvAB = — 19.1 and a = —1.74 are assumed for the 
UV LF at z ~ 9 for our J-dropout selection (using the 
extrapolations from §5.3). The redshift distributions are 
presented in Figure 7. The mean redshifts expected for 
our z-dropout and J-dropout selections are 7.3 and 9.0, 
respectively. The mean redshifts estimated for both se- 
lections are smaller than what we had considered in pre- 
vious works (e.g., Bouwens et al. 2004c; Bouwens et al. 
2005; Bouwens & Illingworth 2006) because of our adop- 
tion of a slightly more inclusive dropout criterion in this 
work. Note that including the likely evolution across the 
redshift range of these selections would likely lower the 
mean redshifts even more (Muhoz & Loeb 2008). 

4.3. Luminosity / SFR density 

We can use the present LF results to refine our con- 
straints of the rest-frame UV luminosity densities at 
z > 7. Previously, we set constraints on these quantities 
based upon a sample of four z ~ 7 z-dropout candidates 
found over the GOODS fields (Bouwens & Illingworth 
2006) and a sample of three possible J-dropout candi- 
dates over the deep NICMOS parallels to the HUDF 
(Bouwens et al. 2005). 

Here we are able to make modest but significant im- 
provements in our measurement of the UV luminosity 
density at z ~ 7 — 10. We make these estimates based 
upon our best-fit LFs at z ~ 7 and z ~ 9. We shall 
consider integrations to a limiting luminosity of 0.2 J* =3 
because this corresponds to the approximate flux limit 
for our z and J-dropout searches (i.e., -£/i6o,AB ~ 28). 



o 
s 



o 
2 



o 



-1.5 
-2 
-2.5 h 



10 5 3 2 



t (Gyr) 
1 



0.6 



0.4 



> 0.2 14.3 



O 
00 

i 

o 

TO 

09 

— 25.5 a 



25 



6 
z 



10 



s 

13 

O 



Fig. 8. — The present constraints on the UV luminosity density 
at z > 7. At z ~ 7, this constraint is shown as a large solid red 
circle while at z ~ 9, it is shown as a la upper limit [red downward 
arrow). These determinations are integrated to 0.2L*_ 3 to match 
the approximate faint limits on our z > 7 galaxy searches. Also 
shown are the determinations of Schiminovich et al. (2005: open 
black squares), Steidel et al. (1999: green crosses), and Bouwens 
et al. (2007: solid red squares) integrated to the same flux limit. 




Fig. 9. — Estimated star formation rate density as a function of 
redshift (integrated down to 0.2 L*_ 3 as in Figure 8). The lower 
set of points give the SFR, density without a correction for dust 
extinction, and the upper set of points give the SFR density with 
such a correction. This is also indicated with the shaded blue and 
red regions, respectively, where the width of these regions show the 
approximate uncertainties estimated by Schiminovich et al. (2005). 
At lower redshift (z < 3), we adopt the dust correction suggested by 
Schiminovich et al. (2005). At z > 6, we adopt the dust correction 
obtained by Bouwens et al. (2006: see also Stark et al. 2007a and 
Stanway et al. 2005) at z ~ 6 from the [/ V-continuum of i-dropouts 
and the Meurer et al. (1999) IRX-/3 prescription. At z ~ 4 and 
z ~ 5, we interpolate between the estimated dust extinctions at 
z ~ 3 and z ~ 6. The symbols are the same as in Figure 8. 

We present these luminosity densities in Figure 8 and Ta- 
ble 6. As is conventional, we also present the equivalent 
star formation rate on this same figure using the Madau 
et al. (1998) and adopting a Salpeter IMF. To show how 
this situation changes when a plausible account of dust is 
included, we plot the equivalent dust-corrected and un- 
corrected SFR densities in Figure 9. We employ the dust 
corrections adopted by Bouwens et al. (2007), which are 
1.4 mag, 1.1 mag, 0.6 mag, and 0.4 mag at z < 3, z ~ 4, 
z ~ 5, and z > 6, respectively. Our use of an evolving 
dust correction from z > 6 to z ~ 3 is motivated by the 
apparent evolution in the t/l^-continuum slope over this 
redshift range (e.g., Stanway et al. 2005; Bouwens et al. 
2006) and the correlation of [/^-continuum slope with 
dust extinction (e.g., Meurer et al. 1999; Reddy et al. 
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TABLE 6 

UV Luminosity Densities and Star. Formation Rate 
Densities. a 







iogi £ 


log 1() SFR density 






(ergs s 1 


(Mq Mpc' 


-3 yr-l) 


Sample 


< z > 


Hz" 1 Mpc- 3 ) 


Uncorrected 


Corrected 15 


z 


7.3 


25.32±0.21 


-2.58 ±0.21 


-2.40 ±0.21 


J 


9.0 


< 25.14 c 


< -2.76 c 


< -2.58 c 


B 


3.8 


26.18±0.05 


-1.72±0.05 


-1.29 ± 0.05 


V 


5.0 


25.85±0.06 


-2.05±0.06 


-1.75 ±0.06 


i 


5.9 


25.72±0.08 


-2.18±0.08 


-2.00 ±0.08 



a Intcgratcd down to 0.2 L* =3 . Based upon LF parameters in 
Table 4 (see §4.3). 

b The adopted dust-extinction corrections are 0.4 mag, 0.6 mag, 
and 1.1 mag at z > 6, z ~ 5, and z ~ 4, respectively. Our use of an 
evolving dust correction from z > 6 to z ~ 3 is motivated by the 
apparent evolution in l/V-continuum slope over this redshift range 
(e.g., Stanway ct al. 2005; Bouwens et al. 2006) and the correlation 
of (7V-continuum slope with dust extinction (e.g., Meurer et al. 
1999; Reddy et al. 2006). 

c Upper limits here are Icr (68% confidence). 

2006). 

5. DISCUSSION 

5.1. How does the volume density of UV -bright galaxies 

evolve? 

Over the past few years, there has been some contro- 
versy regarding how the UV LF evolves at high redshift. 
While some studies have argued that the evolution pri- 
marily occurs at the bright end (i.e., Dickinson et al. 
2004; Shimasaku et al. 2005; Ouchi et al. 2004a; Bouwens 
et al. 2006; Yoshida et al. 2006; Bouwens et al. 2007), 
there have been other efforts which have argued that the 
evolution occurs primarily at the faint end (i.e., Iwata et 
al. 2003; Sawicki & Thompson 2006; Iwata et al. 2007). 
In this work, we found additional evidence to support the 
fact that the most rapid evolution occurs at the bright 
end of the UV LF and that this evolution can be ap- 
proximately described by a change in the characteristic 
luminosity M^y of the UV LF. 

A significant part of the debate regarding the form the 
evolution of the UV LF takes at high redshift is cen- 
tered on what happens at the bright end of the LF (i.e., 
Muv,ab % —20.5). Does it evolve or not? Iwata et al. 
(2007) and Sawicki & Thompson (2006) argue that there 
is little evidence for substantial evolution from z ~ 5 to 
z ~ 3 from their work. However, a quick analysis of re- 
sults in the literature indicate that such evolution is cer- 
tainly very strong, particularly if the baseline is extended 
out to even higher redshifts (z > 6). At z <~ 6, for exam- 
ple, the bright end of the UV LF (i.e., M UV . AB < -20.5) 
has been reported to be deficient by factors of ~ 6 — 11 
relative to z ~ 3 - 4 (Stanway et al. 2003, 2004; Shi- 
masaku ct al. 2005; Bouwens et al. 2006). At z > 7, this 
deficit is even larger, as we can see by comparing the 
number of UV bright galaxies (i.e., Muv.ab < —20.5) 
in our dropout selections with that expected from z <~ 4 
(Bouwens et al. 2007) assuming no-evolution. The com- 
parisons are presented in Tabic 7. At z ~ 7, the vol- 
ume density of UV bright {Muv,ab < —20.5) galaxies 
appears to be 18tiiX smaller than at z ~ 4 (see also 



TABLE 7 

Comparison of the Observed Numbers of Dropouts 
with No-evolution Extrapolations from z ~ 4 a 



Number Evolutionary 
Sample Observed z ~ 4 Prediction* Factor* 3 



Bright 


{M UV , AB < 


-20.5) 


18±H 


z-dropout 


2 


35.9 


J-dropout 





17.9 


> 16 c 


Faint 


(M UViAB > - 


-20.5) 


4i 3 2 


z-dropout 


6 


24.7 


J-dropout 





8.0 


> 7 C 



Extrapolations assume no evolution in the UV LF from 
2 ~ 4 (Bouwens et al. 2007) and were calculated using Eq. 1. 
The sizes and UV colors of [/V-bright galaxies at z > 7 
assumed for these estimates are the same as those given in 
Appendix B. 

b Ratio of volume density of galaxies at z ~ 4 (Bouwens 
ct al. 2007) in some luminosity range to that in our higher 
redshift dropout samples. 

Tower limits here are Icr (68% confidence). 



Bouwens & Illingworth 2006; Mannucci et al. 2007; Stan- 
way et al. 2008), and at z ~ 9, the volume density of UV 
bright galaxies is at least 16 x smaller (68% confidence). 

The only analysis to find mild evolution in the volume 
density of UV bright galaxies from z ~ 7— 10 to z ~ 3 — 4 
is that of Richard et al. (2006) and involves a search for 
z > 6 galaxies behind lensing clusters. However, as we 
argue in Appendix C, it seems likely that a substantial 
fraction of the [/U-bright galaxies in the Richard et al. 

(2006) z ~ 6 — 10 selection are spurious and therefore the 
bright end of their UV LFs is much too high. The evi- 
dence for this is rather striking: not only are none of the 
z > 6 candidates in the Richard et al. (2006) selection 
detected (<2cr) in significantly deeper (~l-2 mag) NIC- 
MOS±IRAC data (11 of their z > 6 candidates), but also 
their reported prevalence is >10x higher than what we 
measure in searches for similar candidates behind lensing 
clusters (Bouwens et al. 2008; see also Appendix C). 

Considering the substantial evolution found from z > 7 
to z <~ 4 in the volume density of UV bright galaxies, it is 
certainly relevant to note that all three previous studies 
(i.e., Iwata et al. 2003; Sawicki & Thompson 2006; Iwata 
et al. 2007) which argued for less evolution at the bright 
end of the LF than at the faint end did not extend to 
redshifts much higher than 4. In fact, the very mild evo- 
lution (or no-evolution) at the bright end of LF claimed 
by the above studies seems to be highly consistent with 
other analyses in this redshift range (e.g., Steidel et al. 
1999; Reddy et al. 2007) which in general find very little 
evolution from z ~ 4 to z ~ 2. Even for the Iwata et al. 

(2007) LF at z ~ 4.7 (which is fairly similar to the UV 
LFs at z ~ 2 — 4 at bright luminosities), one may not 
expect much evolution at the bright end. After all, the 
mean redshift of their bright selection is z ~ 4.5, which 
is not much higher than z <~ 4. 3 

3 The mean redshift for bright galaxies in the Iwata et al. (2007) 
selection is likely to be somewhat lower in redshift than for the 
fainter ones. This is because (1) the Iwata et al. (2007) selection 
window extends to somewhat lower redshift for the redder star- 
forming galaxy population than it does for the bluer star-forming 
galaxy population and (2) brighter star-forming galaxies at z ~ 
4 — 6 are observed to have much redder UV colors in general than 
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TABLE 8 

Best-fit schechter parameters for the UV LFs at 

2 ~ 7 AND Z ~ 9 ASSUMING DIFFERENT MEAN 
C/V-CONTINUUM SLOPES /3'S TO COMPUTE THE SELECTION 
VOLUMES FOR OUR DROPOUT SAMPLES (SEE §5.2). a 



Dropout 


Mean 




<t>* (10~ 3 


Sample 


< z > /3 


IV1 UV 


Mpc~ 3 ) 



z 


7.3 


-2.0 


-19.8 ±0.4 


1 1+ 16 

i - 1 -0.7 


z 


7.3 


-2.5 


-19.8 ±0.4 


1 q+1.7d 
1 -°-0.8 


z 


7.3 


-1.5 


-19.8 ±0.4 


1 q+1.7d 

1 -°-o.a 


J 


9.0 


-2.0 


> -19.6 C 


(i.ir 

(1.1)- 
(l.l)e 


J 


9.0 


-2.5 


> -19. T 


J 


9.0 


-1.5 


> -19. 5 C 



a All LFs assume a faint-end slope a equal to —1.74. 

b See remarks in Table 4. 

Tower limits here are la (68% confidence). 

d Both bluer and redder CV-continuum slopes can re- 
sult in slightly lower estimated selection volumes for star- 
forming galaxies at z ~ 7. 

°Fixed in our fits. 



In summary, it would appear that the most reason- 
able conclusions to draw from these studies are that the 
bright end of the UV LF shows substantial evolution at 
very high rcdshifts (z > 4) and then mild evolution at 
somewhat later times (i.e., from z ~ 4 to z ~ 2). 

5.2. Completeness of current samples in UV -continuum 

slope 

The selection criteria we use to select star-forming 
galaxies at z > 7 (and at z < 6) differ somewhat in 
the efficiency with which they select galaxies that have 
red [/^-continuum slopes 0, i.e., > 0.0. This would 
potentially have an effect on our determinations of the 
UV LFs if such red star-forming galaxies existed in large 
numbers. Fortunately, for these dropout selections, it 
appears that they do not. At z ~ 4, for example, 
both Lyman-break and Balmer-break selections almost 
entirely consist of galaxies with /3's bluer than —1.0 (§4.1 
and Figure 8 of Bouwens et al. 2007; Brammer & van 
Dokkum 2007). At z ~ 5 — 6, the galaxy population 
appears to be bluer yet, with mean /3's of ~ —2 (e.g., 
Stanway et al. 2005; Bouwens et al. 2006; Yan et al. 
2005; Lehnert & Bremer 2003; Hathi et al. 2008), and 
there is little evidence from the present selections that 
the galaxy population is any less blue at z > 7. Indeed, 
the observed Jno — #160 colors arc consistent with a 

the fainter star-forming galaxies at these redshifts (e.g., Meurer et 
al. 1999; Iwata et al. 2007; R.J. Bouwens et al. 2008, in prep). 
As a result, the mean redshift for bright galaxies in the Iwata et 
al. (2007) selection will be lower than for the faint galaxies, and 
therefore the bright end of the Iwata et al. (2007) LF should show 
less evolution than the faint end. 

4 For our z-dropout selection, for example, we would only 
be moderately effective at selecting galaxies with C/V-continuum 
slopes /3 redder than fi = (Figure 2). For our J-dropout selec- 
tion, on the other hand, there is no limit on how red a star-forming 
galaxy can be in terms of its fTV-continuum slope, except as set 
by our restrictions on the IRAC and MIPS colors of sources (i.e., 
3.6/6-5.8// < 1, 3.6/4-24// < 2: see §3.2) which require /3 < 0. For 
the Bouwens et al. (2007) z ~ 4 B-dropout, z ~ 5 V-dropout, and 
z ~ 6 i-dropout selections, we are only effective in selecting star- 
forming galaxies, if their t/V-continuum slopes (3 are bluer than 
= 0.5 (see Figure 8 of Bouwens et al. 2007). 



star-forming population with little dust and 10 8 years of 
constant star formation (where ~ —2: sec Figure Al). 

Of course, even if our dropout selections are largely 
complete in their selection of dropouts over a wide range 
in J7V-continuum slope, it is possible that the UV LF we 
derive may depend somewhat on the distribution of UV 
continuum slopes we are using to model the selection vol- 
umes for our high-redshift dropout samples. We would 
not expect the distribution we are assuming to be too 
inaccurate, as we are able to use this distribution to re- 
produce the observed Jno — -Hi60 colors of our z-dropout 
sample (Appendix B: see Figure Al). Nonetheless, it 
is worthwhile examining the sensitivity of our results to 
small changes in the distribution of UV- continuum slopes 
0. Bouwens et al. (2007) examined the effect that dif- 
ferent assumptions about these slopes would have on the 
UV LF at z ~ 4 — 6 and found changes of < 10% in the fit 
parameters for shifts of 0.7 in the mean i/V- continuum 
slope 0. Repeating our determinations of the UV LF 
at z ~ 7 and z ~ 9 assuming a mean J/V-continuum 
slope of —2.5, we find differences of < 20% in the con- 
straints on M* and (p* (Table 8). Similarly small dif- 
ferences were found relative to previous fits assuming a 
mean [/V-continuum slope of —1.5 (Table 8). This 
suggests that our LF fit results are relatively robust with 
respect to uncertainties about the distribution in t/V- 
continuum slopes 0. 

5.3. Evolution of the LF 

In this work, we derived constraints on the UV LF at 
z ~ 7 and z ~ 9 using both a stepwise and Schechter 
parametrization. For our Schechter parametrization, we 
estimated that the value of Mj} y for the UV LF at z ~ 7 
was equal to -19.8 ± 0.4 mag and cj>* = O.OOllt^ooo? 
Mpc~ 3 . The value of M^y we derived at z ~ 7 is 0.4±0.4 
mag fainter than we found at z ~ 6. This is consistent 
with the trend previously discovered over the redshift in- 
terval z ~ 4 — 6 where the value of M^ v was fainter the 
earlier back in time we probed (see also Yoshida et al. 
2006). We have already remarked that the value of </>* is 
quite consistent with the values found at z ~ 4, z ~ 5, 
and z ~ 6 by Bouwens et al. (2007: see also Yoshida et 
al. 2006). We also attempted to constrain the Schechter 
parameters at z ~ 9 using the results of our J-dropout 
search. Assuming that we can model the changes in the 
UV LF from z ~ 9 using a Schechter function and that 
these changes can parametrized by an evolution in M^ v 
(with <f>* and a fixed, as is suggested by our constraints 
on the UV LF from z ~ 4 to z ~ 7), we found that 
2=10 > —19-6 mag (68% confidence). Plotting the 
constraints on My V versus redshift (Figure 10), it seems 
clear that the observations are consistent with a sub- 
stantial brightening of M^y at early cosmic times. The 
brightening of M^ v is > 1.4 mag from z ~ 9 to z ~ 4. 
This, of course, is in contrast to the behavior of M^ v 
at z < 4 where M^j v reaches a maximum brightness at 
z ~ 2 — 4 and then fades with cosmic time (e.g., Gabasch 
et al. 2004; Arnouts et al. 2005). 

We can put together the constraints on the UV LF 
from z ^ 7 — 10toz^4 (i.e., Figure 5 from this work 
and Figure 3 from Bouwens et al. 2007) to derive approx- 
imate expressions for M- { 



and a as a function of 
redshift (at z > 4). Assuming that the dependence of 
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Fig. 10. — Determinations of the characteristic luminosity M* in 
the rest-frame UV as a function of redshift. Shown are our present 
determination at z ~ 7 (large red circle) and la upper limit at 
z ~ 9 (red downward arrow) . Also included are the determinations 
at z ~ 4, z ~ 5, and 2 ~ 6 from Bouwens et al. (2007: red circles), 
at z ~ 3 from Steidel et al. (1999: green cross), and at 2 < 3 from 
Arnouts et al. (2005: blue crosses), and at z ~ from Wyder et 
al. (2005: blue square). The expected evolution in M* predicted 
from the momentum-driven wind "vzw" model of Oppenheimer 
& Dave (2006) is plotted as the dashed magenta line while the 
evolution expected as a result of a build-up in the mass function is 
shown. The solid black line shows the evolution predicted from the 
halo mass function (Sheth &; Tormen 1999) assuming a mass-to- 
light ratio which is constant while the dotted black line shows this 
evolution assuming a mass- to-light ratio which varies as (1 + 2) — 1 
(which is somewhat less than the (l + z)~ 1,5 scaling of the mass-to- 
light ratio used in Wyithe & Loeb 2006 or Stark et al. 2007c: see 
§5.4). The predicted evolution in M^y is derived from the models 
by examining the changes in the LF (or mass function) at a fixed 
volume density ~ 10 -2 ' 5 Mpc -3 . The above predictions differ 
from those given in Bouwens et al. (2007) due to our use of a more 
current input matter power spectrum (Spergel et al. 2007). The 
observed evolution in M^ v appears to be in reasonable agreement 
with what one might expect from the evolution of the halo mass 
function, assuming an evolution in the mass-to-light ratio for halos 
of ~ (1 + z)- 1 . 

these parameters on redshift can be parametrized in a 
linear way and moreover that the the UV LF can be rep- 
resented by a Schechter function at z > 7 (see discussion 
in §5.5), we derive the following relationships for M£y, 
(/)*, and a for z ~ 4 - 7: 

Mlr V = (-21.02 ± 0.09) + (0.36 ± 0.08)0 ~ 3 - 8 ) 

4>* = (1.16 ± 0.20)10(°- 024±0 - 065 ^- 3 - 8 >10- 3 Mpc- 3 

a = (-1.74 ± 0.05) + (0.04 ± 0.05)0 - 3.8) 

These equations provide us with a convenient way of ex- 
trapolating our observational results back to early times 
(i.e., z > 7). They also provide us with a very quantita- 
tive way of assessing the significance of the evolutionary 
trends identified in My V , (jf , and a. In particular, these 
equations suggest that the observed evolution in Mfj v 
is significant at the 4.5cr level over the interval z ~ 9 
to z ~ 4, while the evolution in cf>* and a is still not 
significant at all (i.e., < la). 

5.4. Relation to the Halo Mass Function 

The rapid brightening of galaxies within the first two 
billion years of the universe is not particularly surprising 
since it is during this time period that L* galaxy ha- 
los (typically 10 10 to 10 12 M Q : e.g., Cooray 2005; Lee 
et al. 2006; Ouchi et al. 2004b; Cooray & Ouchi 2006) 
are largely assembled. The volume density of ~ 10 11 



Mq halos is expected to increase 1000 fold from z ~ 10 
to z ~ 3 — 4 (see, e.g., Figure 2 of Springel et al. 2005). 
Larger increases are expected at > 10 11 Mq, and smaller 
increases are expected at < 10 11 Mq. Thus, if the effi- 
ciency of star formation in galaxies does not decrease 
substantially as a function of cosmic time, we would ex- 
pect substantial increases in the volume density of lu- 
minous star-forming galaxies from z ~ 9 to z ~ 4. We 
would also expect the maximum typical luminosity of 
star- forming galaxies (and thus M^y) to increase over 
this time period (due to the greater evolution expected 
for more massive galaxies) . This is exactly what our LFs 
suggest (Figure 10). 

In light of these expectations, it is relevant to try to 
connect the evolution we observe in the rest-frame UV 
LF to the anticipated changes in the halo mass function. 
Not only does this provide us with a simple test of our 
theory for structure formation, but it also allows us to 
place modest constraints on the efficiency with which ha- 
los in the early universe produce UV photons. 

In our previous work (§5.2 of Bouwens et al. 2007), we 
decided to parametrize the evolution in the UV LF and 
that expected for the halo mass function in terms of a 
single variable for simplicity. Motivated by our obser- 
vational findings (Bouwens et al. 2006; Bouwens et al. 
2007), we took this parameter to be the characteristic 
luminosity for our UV LF and some characteristic mass 
for the halo mass function. The characteristic mass was 
taken to equal the mass where the mass function had a 
comoving volume density of 10~ 25 Mpc -3 to (approx- 
imately) match the volume density of <fi* galaxies. We 
ignored questions about a possible evolution in the shape 
of the UV LF or mass function. 5 Making different as- 
sumptions about how the mass-to-light ratio (halo mass 
to UV light) of halos vary as a function of redshift, we 
developed predictions for how M^ v might be expected 
to change with redshift - assuming a matter power spec- 
trum with erg = 0.9, spectral index n s = 1, CIm = 0.3, 
A = 0.7, and H = 70km/s/Mpc. We found that the 
observations seemed to prefer a scenario where there was 
little evolution in the mass-to-light ratio of halos at early 
times. 

Given the similarity of the present observational results 
to those previously considered by Bouwens et al. (2007), 
one might suspect that the above conclusions would re- 
main unchanged. However, after exploring the sensitiv- 
ity of these conclusions to different inputs, we found that 
they could depend significantly upon assumptions made 
in calculating the halo mass function. In particular, we 
found that for the matter power spectrum preferred by 
the WMAP three year data and the SDSS large-scale 
structure constraints (os = 0.772, n s — 0.948), the evo- 
lution we predicted for the characteristic mass was ~40% 
larger than what we predicted in our previous work using 
er 8 = 0.9 and n s = 1. This translated into a ^40% larger 
predicted brightening of Mfj V with cosmic time than in 
the previous predictions, and so as a result our models 
with no evolution in the mass-to-light ratio (halo mass 
to UV light) were no longer in good agreement with the 

5 We note that we find very similar changes in the so-called 
characteristic mass here if we make the comparisons in terms of 
the average (or even median) mass of halos at the high-mass end 
of the mass function (< 10 -25 Mpc -3 ). 
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observations (see solid line in Figure 10). 

We therefore explored different scalings of the mass- 
to-light ratio of halos to try to regain good agreement 
with the observations. We found reasonable agree- 
ment, adopting a halo mass-to-light ratio that varied as 
(1 + (see dotted line in Figure 10). Such a scal- 

ing would mean that star formation at early times is 
much more efficient than at later times and, in fact, does 
not scale that differently from what one would expect 
based upon the dynamical time scales (which evolves as 
- (1 + z)- 1 - 5 : e.g., Wyithe & Loeb 2006; Stark et al. 
2007c). Of course, it seems prudent to remember that 
these conclusions are likely highly uncertain. They are 
clearly somewhat sensitive to the input matter power 
spectrum that one adopts, and it seems likely that there 
are other uncertainties that will have an effect (e.g., the 
issues discussed in §5.5). 

5.5. Shape of the UV LFs at z > 4 

Of course, the UV LFs we have derived at z > 4 con- 
tain a lot more information than is present in a single 
parameter. They also contain substantial information on 
the overall shape of the UV LF. Perhaps the most no- 
table feature in the shape of the observed LFs at z ~ 4, 
z ~ 5, and z ~ 6 is the presence of a distinct "knee," 
where the LF transitions from a power-law-like behavior 
towards fainter magnitudes to an exponential cut-off at 
the bright end. At lower redshifts, this exponential cut- 
off at the bright end of the LFs is well-known and has 
been hypothesized to occur as a result of a variety of dif- 
ferent factors. Classically, this cut-off has been explained 
as a result of the inefficiency with which gas can cool in 
halos above 10 12 Mq (e.g., Binney 1977; Rees & Ostriker 
1977; Silk 1977) and settle out of hydrostatic equilibrium. 
This explanation has recently been updated and recast 
in terms of a discussion of hot and cold flows (i.e., Birn- 
boim & Dekel 2003; Dekel & Birnboim 2004). The expo- 
nential cut-off at the bright end of the LF has also been 
explained as a result of AGN feedback (e.g., Binney 2004; 
Scannapieco & Oh 2004; Croton et al. 2006; Granato et 
al. 2004), which above a certain mass scale is thought to 
cut off the flow of cold gas to a galaxy and therefore stop 
star formation. One final factor which may be important 
for imparting the sharp cut-off at the bright end of the 
LF, particularly at UV wavelengths, is dust obscuration, 
which becomes substantial (factors > 10) for the most 
luminous systems (e.g., Wang & Heckman 1996). As a 
result, the UV luminosity of star- forming galaxies does 
not tend to reach above some maximum value (see dis- 
cussion in e.g., Adelberger & Steidel 2000; Martin el al. 
2005). 

At very high redshifts, where the typical halo masses 
are smaller, none of the above mechanisms would seem 
to work very well. The typical masses are below those 
where AGN feedback might be expected to be important 
(e.g., Scannapieco & Oh 2004; Dekel & Birnboim 2004) 
or where gas cooling is inefficient (i.e. < 10 12 M Q : Dekel 
& Birnboim 2004). 6 The expected star-formation rates 
at these redshifts would also tend to be much lower and 
therefore dust obscuration would not likely be such an 

6 Nonetheless, it has been argued (e.g., Dekel & Birnboim 2004) 
that at z > 2 gas cooling will occur very efficiently in some regions 
of high mass (> 10 12 Mq) halos through "cold flows." 
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Fig. 11. — Rest-frame UV LFs predicted at z ~ 4, z ~ 5, z ~ 6, 
z ~ 7, and z ~ 9 (dashed blue, green, red, magenta, and blacks 
lines, respectively) using a CLF formalism (see §5.5) compared 
with those from our observations. The stepwise LFs found in the 
observations are presented in the same way as in Figure 4. The 
CLF kernel used to convert mass into light was chosen to fit the 
LF at 2 ~ 4. The predicted LFs using a kernel with no cut-off 
at the bright end of the mass function is shown using the dotted 
line (i.e., assuming L C (M) = (2.51 X 10 22 WH z - 1 )(M/m c ) 1 ' 24 ). 
At z > 4, we adopt the same CLF kernels as at z ~ 4 to calculate 
the UV LF from the halo mass function, but modified to include 
the ~ (l + z)" 1 evolution in the mass-to-light ratio inferred by 
comparing the UV LFs with the expected halo mass function (see 
Figure 10). The LFs predicted using these assumptions do not 
show an extremely prominent "knee" so characteristic of "typical" 
Schechter functions at z > 6. 

important factor either (given the correlation between 
these two quantities: e.g., Wang & Heckman 1996). As 
such, it would seem that there are no obvious mecha- 
nisms to impose an abrupt cut-off at the bright end of the 
UV LF and therefore we might expect that the UV LFs 
at very early times [z > 6) to lack the "knee" so charac- 
teristic of the Schechter parametrization and look much 
more like the mass function. In fact, the UV LFs at z > 4 
derived from a number of recent theoretical simulations 
(e.g., Figure 7 of Finlator et al. 2006; Figure 10 of Oesch 
et al. 2007) have exactly those characteristics. Those 
predicted LFs bear a much closer resemblance to power 
laws than they do to the traditional Schechter function. 
The observation of a significant knee-or cut-off-at the 
bright end of the observed LFs out to z ~ 6 may mean 
that there is some additional physics relevant to the for- 
mation of galaxies that is not fully understood. 

Conditional Luminosity Functions: To illustrate this 
issue quantitatively, we adopt a simplified version of the 
conditional luminosity function (CLF) formalism laid out 
in Cooray & Ouchi (2006: see also Yang et al. 2003; 
Cooray & Milosavljevic 2005). The CLF formalism pro- 
vides us with a way of converting some halo mass func- 
tion to a LF. Typically there is a kernel which when 
convolved with the mass function gives the LF, i.e., 

H L)=fj(L\M)^ d M 
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where 

PYn { log 10 [£/£c(M)] 2 1 
ex f \ 2^ / 

where ^jy is the Sheth-Tormen (1999) mass function, 
where log e 10 « 2.303, where L C (M) gives the t/V lumi- 
nosity of the central galaxy in some halo of mass M and 
where the parameter a expresses the dispersion in the 
relationship between the halo mass and the UV light of 
the central galaxy. Note that we have ignored the contri- 
bution from satellite galaxies to the luminosity function 
in the above equation since they appear to constitute 
< 10% of the galaxies over a wide-range in luminosity 
(see, e.g., Cooray & Ouchi 2006). One aspect about the 
CLF formalism is that it is largely independent of much 
of the normal phenomenology associated with galaxy for- 
mation and therefore it is very clean. There is no attempt 
to explicitly model the large number of physical processes 
likely to be important in the formation and evolution of 
galaxies (which remain a challenge for most semi-analytic 
prescriptions, e.g., Somerville & Primack 1999; Cole et 
al. 2000). 

Making use of the above formalism, we experimented 
with a variety of different paramcterizations to find a 
function L C (M) and parameters a and m c which repro- 
duced the Bouwens et al. (2007) UV LF at z ~ 4. We 

found that L C (M) = (2.51 x 10 22 WHz^j ggg , 

a- = 0.16, and m c = 1.2 x 10 12 M Q worked rea- 
sonably well (see the blue dashed line in Figure 11). 
2.51 x 10 22 WHz -1 is equivalent to -21.91 AB mag. We 
adopted the best-fit values of u 8 , n s , fi m , ^6, an d 
h to the WMAP three year + SDSS data (Spergel et al. 
2007). Our choice of m c set the mass scale where the ex- 
ponential cut-off in the luminosity function takes effect. 
At the high mass end of the mass function, our equation 
for L c asymptotes towards a (M/m c ) ' 2 behavior. The 
physical motivation for this cut-off at the bright-end of 
the mass function would presumably be the same as at 
low redshift, where star formation is not as efficient for 
the most massive halos because of inefficient gas cool- 
ing or AGN feedback. To show the effect the high mass 
cutoff has on the model LF, we also show this LF with- 
out such a cut-off (see the dotted lines in Figure 11). 
At the low-mass end, our equation for L c aymptotes to- 
wards a (M/m c ) 124 behavior. This allowed us to ac- 
count for the fact that the faint-end slope of the UV LF 
is slightly shallower than that of the mass function (i.e., 
— 1.7 vs. —2.0). The chosen parameters yield similar, 
albeit slightly lower, halo masses to those inferred from 
clustering analyses of star-forming galaxies at z ~ 4 (e.g., 
Ouchi et al. 2004b; Lee et al. 2006). 

The question now is whether the above relationship be- 
tween UV light and halo mass (i.e., the kernel (f>{L\M)) 
can reproduce the rest-frame UV LF at z > 4. Given 
the reasonable agreement between the observed evolu- 
tion in My V and that predicted from the halo mass 
function assuming ~ (1 + z) _1 in the mass-to- light ra- 
tio, one might expect that the above kernel may work 
if we multiply the L C (M) expression above to include 
the factor ( pp^g ) ■ With this change, L C {M) equals 

(2.51 x 10 22 WHz- 1 ) fffeQ (^g). ThcUVLFswe 



predict atz~5,z^6,z^7, and z ~ 9 are presented in 
Figure 11 as the green, red, magenta, and black dashed 
lines, respectively. These LFs appear to differ somewhat 
in overall shape from what is found in our observations, 
predicting somewhat fewer galaxies at fainter luminosi- 
ties than are observed, particularly near the "knee" of 
the LF (i.e., at ~ —19.5 AB mag). Most notably, the 
sharp jump in the volume densities of sources observed 
between the bright end of the LF (i.e., ~ —22 AB mag) 
and the "knee" of the LF (i.e., ~ —20 AB mag) is not 
well matched by the predicted LFs at z > 6. 

This effectively illustrates what we were describing 
above. To obtain a LF with a prominent "knee," we re- 
quire some physical mechanism which will truncate the 
LF above some luminosity (or mass) . Since this "knee" of 
the LF would seem to shift to much fainter luminosities 
at early cosmic times (assuming that the evolution of the 
UV LF can be effectively parametrized in terms of My V 
to z ~ 9 as suggested by Figure 10 or Eq. 3), this mech- 
anism needs to become important at progressively lower 
luminosities (and also lower masses), as we move back 
in cosmic time. This is in contrast to the assumptions 
we use in the formalism above, where we have assumed 
that the mass cut-off m c is independent of cosmic time. 
It is not clear that any of the mechanisms proposed to 
produce a cut-off at the bright end of the LF (e.g., AGN 
feedback, dust extinction, inefficient gas cooling at high 
masses) will extend to very low luminosities (or masses) . 
The observations may therefore present us with a bit of 
puzzle. Alternatively, it is also possible given the ob- 
servational uncertainties that the UV LFs at z > 7 will 
increasingly look like the halo mass function and less like 
a Schechter function with a prominent "knee." 

6. SUMMARY 

We have taken advantage of all the deep near-IR and 
optical data that are currently public over and around 
the two GOODS fields (including the HUDF) to search 
for z > 7 z and J-dropout candidates. In total, we make 
use of ~ 23 arcmin 2 of deep NICMOS data and ~ 248 
arcmin 2 of deep ACS+ISAAC/MOIRCS data. Over 
these areas, we identify 8 z-dropout candidates, but no 
plausible J-dropout candidates. By exploring all reason- 
able sources of contamination (T-dwarfs, high-rcdshift 
supernovae, photometric scatter, spurious sources) for 
our modest z-dropout sample, we argue that the ma- 
jority of sources in our selection (> 88%) are bona-fide 
star-forming galaxies at z <~ 7. Using the size-redshift 
trends, UV colors, and spatial profiles of star-forming 
galaxies at z ~ 4 — 6, we make plausible estimates of 
the selection volume for identifying z ~ 7 z-dropouts 
and z ~ 9 J-dropouts. We then use these samples and 
selection volumes to derive the LF at z <~ 7 and set con- 
straints on the LF at z ~ 9. We consider both stepwise 
and Schechter paramcterizations for the UV LF. 

Our primary conclusions are: 

(1 ) UV LF at z ~7: The rest-frame UV LF shows sig- 
nificant evolution from z <~ 7 to z ~ 6 (99% confidence). 
Our best-fit <j>* and M^y for the UV LF at z ~ 7 is 
0.0011^;^ Mpc~ 3 and -19.8 ± 0.4 mag, respectively, 
adopting the same value for the faint-end slope a (—1.74) 
that we previously found at z ~ 6. This value of M^ v is 
nearly 1.2 ± 0.4 mag fainter than the value of M^ v that 
has previously been found at z ~ 4 and 0.4 ± 0.4 mag 
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fainter than the value found at z ~ 6. By contrast, the 
value of </>* that we determined at z ~ 7 is not signifi- 
cantly different from those values determined at slightly 
later times (z ~ 4 — 6: Bouwens et al. 2007). Similar to 
several previous studies (Bouwens et al. 2006; Yoshida et 
al. 2006; Bouwens et al. 2007), this suggests that the evo- 
lution of the UV LF at z > 4 can largely be described by 
an evolution in the characteristic luminosity Mjtjy, with 
no significant changes required in the other variables. 

(2) Constraints on the UV LF at z ~ 9: Since we did 
not find any credible J-dropout candidates in our search 
fields, we were not able to derive best-fit parameters for 
the UV LF at z <~ 9 and could only set constraints on 
<j>* and M^y. To interpret the implications of our J- 
dropout search results for the UV LF at z ~ 9, we as- 
sumed that the value of <f>* and a there were roughly 
equal to the values (i.e., 0.0011 Mpc -3 and —1.74) found 
at z ~ 4 — 8. We then derived lower limits on the value 
of My V by using our J-dropout search results and the 
effective selection volumes we calculated for z ~ 9 galax- 
ies using numerous Monte-Carlo simulations. We found 
that My V > —19.6 mag at 68% confidence. These search 
constraints are inconsistent with no evolution from z ~ 9 
to z <~ 7 at 80% confidence. 

(3) Evolution in the volume density of UV -bright 
galaxies: Comparing the no-evolution predictions of the 
Bouwens et al. (2007) z ~ 4 UV LF at bright magni- 
tudes with those observed in the ISAAC, MOIRCS, or 
NICMOS data (see Table 7: see §5.1), we estimated that 
the bright end {M UV . AB < -20.5) of the UV LF is at 
least 18^10 x lower at z ~ 7 than it is at z ~ 4. Making 
a similar comparison between the number of J-dropouts 
predicted assuming no evolution from z ~ 4 and the zero 
candidates found in our NICMOS search fields, we esti- 
mated that the bright end of the LF (M uv ,ab < -20.5) 
is at least 16 x lower at z ~ 9 than it is at z ~ 4 (68% 
confidence). This provides strong evidence against evo- 
lutionary models with little evolution at the bright end 
of the UV LF. 

(4) Empirical fit to the UV LF from z <~ 9 to z ~ 4: 
Putting together our constraints on the UV LF from 
z ~ 9 to z ~ 4, we derived fitting formula to describe the 
evolution of Mfry, <fi*, and a as a function of redshift. 
The results are M* y = (-21.02±0.09) + (0.36±0.08)(z- 
3.8), <j>* = (1.16±0.20)10(°- 024±0 - 065 )( z - 3 - 8 )lO- 3 M P c- 3 , 
and a = (-1.74 ± 0.05) + (0.04 ± 0.05) (z - 3.8). These 
formulae imply that the evolution in Mj!j v is significant 
at the 4.5ct level. The evolution in 4>* and a is not yet 
significant. 

(5) Relationship between the evolution of the UV LF 
at z > 4 and that expected for the halo mass function: 
The evolution we observe in the characteristic luminos- 
ity My V versus redshift matches the evolution expected 
for the halo mass function, assuming that the mass-to- 
light ratio (halo mass to UV light) of halos evolves as 
- (1 + z)- 1 with cosmic time. This finding is somewhat 
different from what Bouwens et al. (2007) inferred about 
the mass-to-light ratios of halos based upon very similar 
observational constraints. The reason for the change is 
our use of more current constraints on the matter power 
spectrum (i.e., erg = 0.772, spectral index n s = 0.948: 
Spergel et al. 2007) rather than the a$ = 0.9, n s = 1 
power spectrum we were assuming in Bouwens et al. 



(2007). Effectively, this illustrates how sensitive galaxy 
evolution can be at high-redshift to the input matter 
power spectrum. 

(6) Shape of the UV LF at z > 4: At redshifts as early 
as z ~ 6, the UV LF is remarkably well-represented by a 
Schcchter function (Bouwens et al. 2007). The position 
of the "knee" of the LF is known to < 0.2 mag at z < 6 
(Bouwens et al. 2007) and would seem to be known to 
< 0.4 mag at z ~ 7. However, it is unclear why the UV 
LF would show a prominent knee at such early times. 
Most of the physical mechanisms thought to create this 
knee in the LF at lower redshifts (e.g., AGN feedback or 
a post-virialization cooling condition) would not seem to 
be effective at the lower masses relevant at these times. 

The constraints we have on the UV LF at z > 7 are al- 
ready quite notable, especially given the significant chal- 
lenges that exist to the detection of z > 7 galaxies with 
current observational facilities. Putting together all the 
data now available, we can plausibly determine the vol- 
ume density of £/V-bright galaxies at z ~ 7 to within a 
factor of 2 and set meaningful constraints on the volume 
density of similarly luminous galaxies at z ~ 9. Com- 
paring the present search results with those obtained at 
lower redshift (z ~ 4 — 6), we observe a rather substan- 
tial increase in the volume density of most luminous UV- 
bright galaxies, from redshifts z>7toz~4 - suggest- 
ing that such luminous and likely prodigious star-forming 
systems only start to become common relatively late in 
cosmic time (i.e., 1.5 Gyr after the Big Bang). These 
galaxies are presumably built up through the merging 
and coalescence of lower mass galaxies. 

Over the next year, we expect our constraints on the 
bright end of UV LF at z > 7 to improve substantially 
as a result of both ground and space-based imaging pro- 
grams. From the ground, these improvements will come 
as a result of the deep wide-area imaging being obtained 
over the GOODS fields with wide-area near-IR imagers 
such as MOIRCS on Subaru telescope and HAWK-I on 
VLT-Yepes. From space, these improvements will come 
from two large-to-medium-sizcd HST NICMOS programs 
(GO-11082 and GO-11144 are surveying ~ 60 arcmin 2 to 
a depth of ~ 26.5 AB). Of course, for truly significant 
advances in the size and depth of current samples, we will 
need to wait until late 2008 when WFC3 with its pow- 
erful IR channel will be installed on the Hubble Space 
Telescope. 
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Fig. Al. — Comparison between the observed properties of the i and z dropouts (black histogram) and that expected based on our 
model for star- forming galaxies at z > 6 (red lines: see §4 and Appendix B for a description of this model). The leftmost panels show 
comparisons between the observed and expected half-light radii (r^;) while the rightmost panel shows comparisons between the observed 
and expected UV colors. Only galaxies brightward of -Hieo.AB = 26.4 and -Hi60,AB = 27.0 in our HUDF and NICPAR selections are 
included in the distribution of half-light radii. This leaves us with some leverage in examining the size (surface brightness) distribution. 
Otherwise, selection effects would tend to dominate, excluding all the larger sources from these comparisons and we would be left with 
only those sources whose sizes were very similar to that of the PSF. For the comparisons we make involving UV colors, all search fields are 
considered. In general, good agreement is found between the observed and expected properties, providing us with a reasonable amount of 
confidence in the model we are using to calculate the effective selection volumes for star-forming galaxies at z > 7. 
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APPENDIX 

A. LARGE SCALE STRUCTURE UNCERTAINTIES 

In this work, we derived our best-fit LF at z ~ 7 by comparing the observed surface density of galaxies with what 
we would expect based on a model LF and then maximizing the likelihood of the model LF. We assumed Poissonian 
statistics in computing the likelihood of various model LFs (§4). Our decision to adopt a Poissonian model in computing 
these likelihoods rather utilizing a more conventional approach like SWML was to take advantage of the substantial 
amount of deep, wide-area near-IR imaging (e.g., the ISAAC or MOIRCS data) we have to constrain the UV LF 
at z ~ 7. Since most of these imaging data contain no z ~ 7 candidates, we would not be able to use that data 
for improving our determination of the UV LF at z ~ 7, had we chosen to use the conventional SWML or STY79 
techniques. However, modelling our observational results using Poissonian statistics (and assuming no large-scale 
structure variations), we are able to take advantage of this information to constrain the shape of the UV LF. 

Of course, in modelling the observational results using this approach, we are inherently assuming that the volume 
density of high-redshift star- forming galaxies is statistically homogeneous across the sky. In reality, however, the 
distribution of galaxies is modestly clustered, and therefore the mean volume density of galaxies will show significant 
variations as a function of position. These variations in the mean density occur as a result of large-scale structure and 
will have an effect on the best-fit UV LF we derive at z ~ 7, scattering it away from the cosmic average. 

To estimate the effect that such variations will have on our fit results, we ran a set of Monte-Carlo simulations. We 
varied the surface density of model dropouts in each of our search fields by ~ 30% RMS (see §4.1) and then recomputed 
the UV LF in the same way as in §4.2. We found la variations of ±0.1 mag in the values of My V and ±30% in the 
values of cf>* at z ~ 7. Adding these uncertainties in quadrature with those deriving from Poissonian uncertainties 
(§4.2), the total uncertainty in <f)* at z ~ 7 is 0.4 dex and My V is ±0.4 mag. 

The effect of large-scale structure variations on the constraints we place on the UV LF at z ~ 9 is much smaller in 
general. Given that we have no sources in our J-dropout sample, the Poissonian uncertainties (± ~ 100%) dominate 
over the expected field-to-field variations, which are ± ~ 20% (see §4.2 or Bouwens et al. 2005). 
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B. PLAUSIBLITY OF OUR MODEL FOR ESTIMATING THE SELECTION VOLUMES OF 

DROPOUT SAMPLES AT Z > 7 

Besides small number statistics, perhaps the biggest uncertainty in the present determination of the rest-frame UV 
LF at z > 7 is in our estimates of the selection volumes for the two dropout samples. We derived these estimates using 
the following model for the sizes, morphologies, and UV colors of galaxies at z > 7. We assumed that star-forming 
galaxies at z > 7 have similar pixel-by-pixel profiles to the z ~ 4 i?-dropouts in the HUDF, but with physical sizes 
scaled by a factor of (1 + z) -1 - 1 (e.g., at z ~ 7 galaxies of a given luminosity have approximately half the physical size 
of star- forming galaxies at z ~ 4 and therefore four times the surface brightness). The [/U-continuum slopes f3 were 
taken to have a mean of —2.0 and la scatter of 0.5. Both assumptions were chosen to match the trends found in the 
observations from z ~ 6 to z ~ 2 for the sizes (e.g., Bouwens et al. 2006; Bouwens et al. 2004a; Ferguson et al. 2004) 
and colors (Stanway et al. 2005; Bouwens et al. 2006). 

Of course, just because the present model is a reasonable extrapolation of the trends (in sizes and colors) found at 
lower redshift docs not mean it will accurately represent star-forming galaxies at z > 7. Therefore, it makes sense 
for us to compare this model with the limited data on z > 6 galaxies where they exist. We begin by comparing the 
observed J — H colors of z-dropouts in our sample with that predicted using our best-fit LF (<f> ~ 0.0011 Mpc~ 3 , 
M^ 9m AB ~ —19.8 mag, a ~ —1.74) and the model we just described above (see also §4) for calculating the selection 
volumes. We include the comparison in the rightmost panel of Figure Al. 

Similarly, we can compare the half-light radii of z-dropouts in our sample with that predicted from our models. 
To ensure that the size distribution we derive for z-dropouts is not dominated by surface brightness selection effects 
(which only allows for the selection of point-like sources near the magnitude limit), we only consider the sizes of 
z-dropouts that are > 0.7 mag brighter than the magnitude limit. Unfortunately, this only leaves us with ~ 2 galaxies 
for making this comparison. The result is shown in the middle panel of Figure Al. The model predictions appear to 
be in reasonable agreement with the observations. As one additional check on the size-redshift scaling we are using to 
model the sizes (surface brightnesses) of z ~ 7 galaxies, we also include a comparison of our model predictions with 
the half-light radii of z ~ 6 i-dropouts in our search fields (leftmost panel of Figure Al). 

C. COMPARISON WITH THE RICHARD ET AL. (2006) RESULTS 

In the present work (e.g., §5.1 or Table 7), we found a significantly smaller (i.e., ~ 1/18) volume density of galaxies 
at the bright end of the UV LF at z > 7 than has been inferred at z ~ 4 (e.g., Steidel et al. 1999; Reddy et al. 2007; 
Bouwens et al. 2007). This is in significant contrast with the volume densities of z > 6 galaxies inferred by Richard et 
al. (2006) from searches for z and J dropouts in ~ 12 arcmin 2 of deep (~ 24.5 AB mag at 5cr) near-IR data behind two 
low-redshift lensing clusters Abell 1835 and AC114. Richard et al. (2006) found 13 first and second category z and J 
dropouts in their search that they consider plausible z > 6 galaxies. Modelling the lensing amplication, contamination, 
and overall incompleteness of their search, Richard et al. (2006) estimated the UV LFs at z ~ 6 — 8 and z ~ 8 — 10. 
In significant contrast to our results, the volume density of bright galaxies in the Richard et al. (2006) LFs were found 
to be very similar to what has been found at z ~ 3 — 4 (i.e., Steidel et al. 1999; Bouwens et al. 2007). 

What is the reason for these substantial differences? While one possible explanation could be "cosmic variance" (i.e., 
field-to- field variations) , simple estimates for the variance one expects in the volume density of high-redshift sources 
behind a typical cluster is of order ~50% RMS, and therefore much smaller than the factor of > 10 differences seen 
here. These "cosmic variance" estimates assume a pencil beam geometry, a redshift window of width Az = 1, an 
angular diameter of ~ 30 arcmin 2 , and a bias of ~ 7, which is appropriate for galaxies at z ~ 7 with a volume density 
of - 10~ 3 - 5 Mpc" 3 (see also Somerville et al. 2004; Trcnti & Stiavelli 2008). 

Instead a much more probable explanation for the discrepancy seen here is simply that a substantial fraction of the 
Richard et al. (2006) z > 6 candidates are spurious in nature. After all, none of their best candidates were detected 
at significance levels much higher than 2.5er, making each of the candidates in their sample open to question. Clearly, 
given the area of empty sky within their 12 arcmin 2 search area, it would not be difficult to imagine a small number of 
these empty areas of sky showing large enough positive fluctuations to satisfy their detection thresholds just by chance! 
To their credit, Richard et al. (2006) recognized the significance of this concern and attempted to model this effect. 
After extensive Monte-Carlo simulations, they argued that the number of dropouts found in their samples significantly 
exceeded that expected from the noise. Unfortunately, such techniques can be highly sensitive to even the smallest 
errors in the noise model, and therefore it is not entirely clear how much weight we should give to the simulations 
Richard et al. (2006) have run to control for this process. 

A much more straightforward and robust way of estimating the spurious fraction is to look at their candidates in 
deeper near-IR data (as Bremer et al. 2004 used to evaluate the reality of the Pcllo et al. 2004 z ~ 10 candidate). 
Such data are available for two deep ~ 0.7 arcmin 2 NIC3 i/i6o exposures over each of their two clusters (AC114 and 
Abelll835). These NICMOS pointings reach to ~ 26.8 AB mag (5 a), ~ 2 mag deeper than the ISAAC data used by 
Richard et al. (2006) in their searches. These pointings provide very deep i?i6o-band coverage for 2 first and second 
category z > 6 candidates reported by Richard et al. (2006) and 5 first, second, and third category candidates. The 
first and second category candidates of Richard et al. (2006) were the z > 6 candidates to which they ascribed the 
highest confidence while the third category candidates were ascribed slightly lower confidence. Of the 5 such candidates 
with very deep i?i6o-band coverage, none show up at > 2a significance, strongly arguing that a substantial fraction of 
their candidates do not correspond to real sources. Taken by itself, it implies that > 80% of the Richard et al. (2006) 
first, second, and third category z > 6 candidates are spurious (68% confidence). If we only consider first and second 
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category candidates (the highest confidence candidates) with deep NICMOS data, the NICMOS non-detections imply 
that > 56% of the Richard et al. (2006) candidates are spurious (68% confidence). 

We also do not find it particularly encouraging that none of the dropout candidates in the Richard et al. (2006) 
sample are detected in the deep IRAC 3.6/u imaging which exist around these clusters (Schaerer et al. 2007). There are 
six sources in the Richard et al. (2006) sample which arc sufficiently isolated from the neighbors for these upper limits 
to have particular significance. Given the reported _ff-band magnitudes of the Richard et al. (2006) z > 7 candidates 
(ranging from 24.6 to 25.2 AB mag), we would suspect it should be easy to detect these sources in the very deep 
IRAC data available (reaching to 25.8 AB mag at 2a), especially given that typical (zsso — 3. 6{im)AB colors for other 
z ~ 6 — 8 candidates range from —0.5 mag to 1.5 mag (Yan et al. 2006; Labbe et al. 2006; Eyles et al. 2007; Bradley 
et al. 2008). While Schaerer et al. (2007) suggest that the IRAC non-detections argue for a very young population, 
we do not think the candidates' reality is well enough established to draw any conclusions about the properties of the 
purported sources. Rather we find it worrying that none of the available NICMOS or IRAC data over these clusters 
have provided us with evidence that any of the Richard et al. (2006) candidates are real! 

Finally, we note the overall prevalence of z > 6 candidates found by Richard et al. (2006) behind AC114 and Abell 
1835 is at least > 10 x higher than found in searches around other lensing clusters with much deeper, higher quality 
NICMOS data. Bouwens et al. (2008) have quantified this by examining the deep (> 26.5 AB mag at 5<r) J and 
#-band imaging data around 6 galaxy clusters (CL0024, MS1358, Abell 2218, Abell 2219, Abell 2390, Abell 2667) 
and the deep (> 26 AB mag at 5a) J-band imaging around three other galaxy clusters (Abell 1689, Abell 1703, IE 
0657-56). In total, the deep near-IR imaging data around these clusters cover > 20 arcmin 2 in total, almost double 
what was considered in Richard et al. (2006) and extending some > 1 — 2 mag deeper. Only one z > 7 candidate (i.e., 
Al689-zDl: Bradley et al. 2008) and one z <~ 6.5 candidate (Kneib et al. 2004) is found over these fields to the same 
magnitude limit (i.e., 25.5 AB mag) as is considered in the Richard et al. (2006) search. This is to be compared with 
13 first and second category z > 6 candidates found by Richard et al. (2006) in a search over just half the area. Since 
it seems unlikely that high-rcdshift dropouts are > 10 x more common behind the two particular clusters considered 
by Richard et al. (2006) than the 9 other clusters with deep NICMOS data, this again suggests that the majority (i.e., 
> 90%) of the z > 6 candidates found in the Richard et al. (2006) search are simply spurious sources. 

D. Z > 7 DROPOUTS IN PREVIOUS SAMPLES NOT IN CURRENT SELECTIONS 

There were a small number of candidates in our previous samples (3 z-dropout candidates from Bouwens et al. 2004c 
and 3 J-dropout candidates from Bouwens et al. 2005) which did not make it into our current selections. In this 
section, we describe why those particular candidates are not included in our current selections. 

Let us first consider the three z-dropout candidates from the HUDF (Bouwens et al. 2004c) which do not make it 
into our current samples. As has already been discussed in Bouwens & Illingworth (2006) and Labbe et al. (2006), two 
of those candidates (UDF-818-886 and UDF-491-880) were the direct result of the "Mr. Staypuft" anomaly (Skinner 
et al. 1998), which produces faint electronic ghosts 128 pixels from bright stars. Since the "Mr. Staypuft" anomaly 
is now carefully modelled and explicitly removed in our "nicred.py" reductions (Magee et al. 2007), this anomaly 
should not pose a problem for our current selections. The third candidate (UDF-825-950) may still be a high-redshift 
star- forming galaxy, but did not match the S/N thresholds used for the current selection. 

Now we consider the three J-dropout candidates reported by Bouwens et al. (2005). What became of these three 
candidates? Two of the three candidates (UDFNICPAR1-04151142 and UDFNICPAR2-09593048) show detections in 
the ultra-deep U-band imaging recently taken over the NICMOS parallels to the HUDF and therefore are clearly not 
z ~ 9 — 10 galaxies. The other candidate (UDFNICPAR1-05761077) is not detected at high enough significance in our 
current reductions of the NICMOS parallels to the HUDF and appears to have been spurious. Spurious sources are a 
particular concern over the NICMOS parallels to the HUDF given the minimal amounts of dithering used in acquiring 
the data. Therefore, substantial care has been taken to accurately characterize the read-out properties of individual 
pixels and correct for any pixel-by-pixel detector (or reduction) signatures (see §2.1). 

The loss of these three z ~ 9 — 10 J-dropout candidates from our samples should not be a significant cause for 
concern. In the Bouwens et al. (2005) analysis, the emphasis was on setting an upper limit on the volume density of 
C/U-bright z ~ 9 — 10 galaxies and therefore we tried to err on the side of including candidates in our samples if there 
was any uncertainty about their nature. 



